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ABSTRACT 

The young stellar cluster illuminating the W40 H II region, one of the nearest 
massive star forming regions, has been observed with the ACIS detector on board 
the Chandra X-ray Observatory. Due to its high obscuration, this is a poorly- 
studied stellar cluster with only a handful of bright stars visible in the optical 
band, including three OB stars identified as primary excitation sources. We 
detect 225 X-ray sources, of which 85% are confidently identified as young stellar 
members of the region. Two potential distances of the cluster, 260 pc and 600 pc, 
are used in the paper. Supposing the X-ray luminosity function to be universal, 
it supports a 600 pc distance as a lower limit for W40 and a total population of at 
least 600 stars down to 0.1 M under the assumption of a coeval population with 
a uniform obscuration. In fact, there is strong spatial variation in i^ s -band-excess 
disk fraction and non-uniform obscuration due to a dust lane that is identified in 
absorption in optical, infrared and X-ray. The dust lane is likely part of a ring of 
material which includes the molecular core within W40. In contrast to the likely 
ongoing star formation in the dust lane, the molecular core is inactive. The 
star cluster has a spherical morphology, an isothermal sphere density profile, and 
mass segregation down to 1.5 M . However, other cluster properties, including 
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a < 1 Myr age estimate and ongoing star formation, indicate that the cluster 
is not dynamically relaxed. X-ray diffuse emission and a powerful flare from a 
young stellar object are also reported. 

Subject headings: ISM: individual objects (W40) - open clusters and associa- 
tions: individual (W40) - protoplanetary disks - stars: pre-main sequence - stars: 
formation - X-rays: stars 



1. Introduction 



The W40 (jWesterhoutl Il958l ) complex is a star-formation regio n in the constellation 



naver 



1970) 



20051 ) with 



Serpens. The complex is made up of the molecular cloud G28.8 +3.5 (IGoss fc S 
with a diameter of ~ 1°, the dense molecular core DoH 279-P7 (IDobashi et al. 
a diameter of ~ 20', an OB star cluster, and an H II blister region. The complex is behind 
a screen of Ay > 5 mag of extinction, so only a handful of stars can be seen in the optical 
band (Figure [lb). Because of the high absorption, the stellar cluster has been little studied. 
The cluster is relatively nearby, but the distance has not been sufficiently constrained and 
estimates vary from 300 pc to 900 pc. Pr evious obser vations show that W40 is partially 
embedded in its natal molecular cloud (e.g. IValledll987l ). which indicates that the cluster is 



no more than a few million years old. 

Infrared (IR) has been the primary wavelength for studies of star formation. However, 
certain abilities of X-ray studies make observations particularly useful for examining young 
stellar clusters (YSCs). X-rays effectively distinguish cluster members from foreground or 
background stars since young stellar o bjects have typical X- ray luminosities 10 3 -10 5 times 
greater than main-sequence stars (e.g.. iPreibisch et al.ll2005l ). The high penetrating power 
of hard X-rays is required to study young stars which are often embedded in na tal molecular 
material. The Chandra Orion Ultradeep Projec t (COUP; iGetman et al.ll2005l ) reveals stars 
with up to 500 mag of visual-band extinction (jGrosso et al. 2005 ). Another advantage of 
studies in the X-ray is that bright nebulosity, which may be a source of confusion in the 
IR or optical, is avoided at X-ray wavelengths, yielding greater source detection sensitivity. 
Differences between X-ray brightness of T-Tauri and OB stars are not so dramatic as in 
optical and in IR bands, so X-rays often allow detection of low mass pre-main sequence (PMS) 
stars close to OB stars. X-ray observations also preferentially detect pre-main-sequence 
(PMS) stars without disks (Class III), unlike the IR, where stars possessing circumstellar 
disks (Class II or 1/0) are preferentially detected due to IR excesses. Thus, X-ray observations 
reduce biases present in IR surveys of young stellar objects (YSOs), and X-ray and IR surveys 
are complementary to each other. 
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The ability to observe stars with large optical extinctions is essential in any study of 
W40 due to the large amount of obscuring material along the line of sight. The high spatial 
resolution of Chandra is also necessary to observe stars in the center of the cluster. The 
X-ray is used to select cluster members for our analysis of W40, and stars with masses down 
to 0.1 M are detected with the data from the 38 ksec Chandra observation. The analysis of 
sources of contamination for X-ray observations of stellar clusters and effective methods to 
pick out cluster stars from samples t hat include foreground stars and extragalactic objects 



are described by lGetman et al.l (120061 ). Once a list of W40 cluster members has been created, 
both X-ray and other wavelengths are used to determine properties of stars in the sample. 
In addition, the spatial distribution of the cluster stars is analyzed. Other interesting phe- 
nomena are present in the X-ray data, including diffuse emission and a large flare from an 
intermediate mass young star. 

The observation and source lists are described in £j2j Optical, infrared, and radio coun- 
terparts are identified in £J3] and near-IR (NIR) stellar properties are derived. Stellar and 
extragalactic contaminants in our source lists are identified in §3.31 X-ray properties are 
derived in §H Confirmatory results relating X-ray and IR properties are presented in §0 
The X-ray luminosity function (XLF), cluster distance, total cluster population, intrinsic 
disk fraction, and age are discussed in £j6l The initial mass function (IMF) is given in £0 
Massive star candidates are discussed in £JBJ Spatial structure is discussed in j|9] in which the 
detection of an absorbing filament, spatial segregation of high versus low- mass stars, and 
spatial inhomogeneity in disk fraction of stellar populations is reported. In §nUl we end with 
the analysis of X-ray diffuse emission and mention a powerful flare from an intermediate 
mass young star. 



1.1. Past Studies of W40 



T he distance to W40 has been controversial. Radial velocity based distance measure- 



ments (IReifenstein et al.lll970l ; |Powneslll970l ; iRadhakrishnan et al.lll972l ) usin g the ISchmidt 



( 1965) model of the galaxy collectively suggest a distance of 600 ± 300 pc ( Vailed 1987 ). 
Shaver fc Gosd (119701 ) also find a distance of 7 00 pc using galac t ic kin ematics. Other dis 



tance measurements fall into this range too. iCrutcher fc Chul (119821) find a dist ance of 
400 pc assuming a spectral type for the brightest star of B0V. ISmith et al.l (119851 ) find a 
distance of 700 pc using measurem ents of W40's stellar component in radio/infrared con- 
tinuum. iKolesnik fc Iurevichl (jl983[ ) use OH absorption-line measurements of the cloud and 
calculate a distance of 600 pc. T hese estimates range f rom 4 00 pc to 700 pc, and we adopt 
a distance of ~600 pc in line with lRodney fc Reipurthl (120081 ) and investigate the validity of 
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this distance in the analysis of the Chandra data. W40 has Galactic coordinates I = 28.8°, 
b = 3.5°. The fairly high Galactic latitude would place the cluster more than 50 pc above 
the Galactic plane if the distance were greater than 800 pc. 



The Aquila Rift, which contains the young star cluster Serpens South (jGutermuth et al. 



20081 ). is projected near W40 in the sky (Serpe ns South is 0.5° to the west; see Figured^,) and 
is only 260±37 pc distant ( jStraizys et al.lll996l ). which is much closer than typical distance es- 
timates for W40. However, W4 and Serpens South are usually regarded a s separate objects 



due to their different distances (ptraizys et al.ll2003 



Gutermuth et al. 



Zhu et al 



200C 



2008h. and 600-700 p 



(]; Zinnecker Sz Yorke 2007 



is the most commonly cited distance f or W40 (cf. 
Rodriguez. Rodney, fc Reipurth 1120101 ) . The nearby Serp ens Main cluster is 3° north of W40 
and is at a third distance of 415 ±15 pc (|Dzib et al.ll2010r ). established through parallax mea- 
surements of its brightest star. iDzib et al.l (120101 ) notes that the superposition of the three 
distinct star-forming regions is not unexpected since the line of sight passes just above the 
Galactic plane. Alternatively, iBontemps et al.l ( 120101 ) assume that W40, the Aquila Rift, and 
Serpens Main are parts of the same star-forming region and are therefore at the distance of 
260 pc for the Aquila Rift. This paper will comment on implications of this lower distance 
for W40. 



The core of the W40 molecular cloud was mapped by lZeilik fc Ladal (I19781) in 12 CO an d 
13 C O emissions. The mas s of th e core is estimated to be ~ 200 — 300 M by lzhu et al. ( 2006 ). 
and iRodney fc Reipurth! (120081 ) report a probable mass for the entire cloud of ~ 10 4 M . 
The center of the star cluster is l ocated at the east edge of the molecular core. Three 



probable OB stars are identified by lSmith et al.l (119851 ) as important ionizing sources for the 
cluster. There is e vide nce for circumstella r mate rial around two of these stars detected by 
Smith et al.l (119851 ) and lVallee fc MacLeod! (119941 ) . An observation of W40 by the Midcourse 
Space Experiment (MSX) reveals the H II region to have an hourglass shape, oriented with 
one lobe to the southeast and another to the northwest, 17' x 30' in extent (see Figure 6 in 
Rodney fc Reipurth! 120081) . Th e clust er core is located just to the northwest of the waist of 



the hourglass. IBontemps et al.l (120101 ). who assume that W40 and Serpens South are part of 
the same star-forming complex at a distance of 260 pc, report ^20 protostellar can didates in 



Herschel observations. A weak outflow is detected in CO lines by IZhu et al.l (120061 ) centered 
on the molecular core, although the driving source is unseen. 
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Chandra Observation and Source List 



2.1. Observation and Data Reduction 



The Chandra X-ray Observatory (CXO) observed W40 for ~38 ks during August, 2007 
using the ACIS camera (IGarmire et al.ll2003l ). The observation has a nominal pointing of 
R.A.= 18 h 31 m 23?8 and decl.= -2° 05' 59" (J2000), the approximate center of the W40 
cluster, and a roll angle of 251 degrees. Data from the imaging array (ACIS-I) are ana- 
lyzed. The ACIS-S3 chip was also operational, but these data are omitted since the point 
spread function is degraded due to the large off-axis angle. The data were taken in Very 
Faint telemetry mode to improve sensitivity to faint sources. No variations in background 
levels are detected. Properties o f the telescope and detector are available via the Proposers' 
Observatory Guide Jcxdl2009h . 



A data reduction protocol that has been used f or Chandra ACI S -I obs ervations of other 
young-star clusters is followed (see Appendix B of iTownsley et al.l (120031 ) for more details 
on the steps used for data reduction). Briefly, using the tool acis-process-events from the 
CIA03.4 software package (IFruscione et al.ll2006l ). the latest calibration information on time- 
dependent gain and a Penn State bad pixel mask are applied, background event candidates 
are identified in the very faint mode, and the data are corrected for CCD charge transfer 
inefficiency Using the acis-detect_afterglow tool, additional afterglow events not detected 
with the standard Chandra X-ray Center (CXC) pipeline are flagged. The event list is 
cleaned with the "grade" (only ASCA grades 0,2,3,4,6 are accepted), "status"§, "good-time 
interval", and energy filters. A small astrometric correction is applied based on matches 
between the eight brightest Chandra sources with small off-axis angles and optical /infrared 



source s from the Naval Observatory Merged Astrometric Dataset (NOMAD; IZacharias et al. 



2004al ). The slight point-spread function (PSF) broadening from the CXC software position 
randomizations was removed. 



2.2. Source Detection and Characterization 

Figure [2b, d shows the ACIS-I field in which more than 200 X-ray sources are seen. 
Candidate X-ray sources in the Chandra data are detected using an aggressive method that 
minimizes missed sources but allows spurious sources. Data images and exposure maps with 



1 As some of the background and afterglow event candidates from bright sources are real they are retained 
for our spectral analysis, and the non-controversial status=0 filter is applied only to build the event list for 
source search. 
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spatial resolutions 0.5" and 1" pixel -1 are used for source detection on the center of ACIS-I 
and the full ACIS-I, respectively. A wavelet-based source detection procedure (wavdetect, 
Freeman et al.l 120021 ) is run with a threshold of P=10 -5 . Sources are detected on a soft 
(0.5-2.0 keV), a hard (2.0-8.0 keV), and the total (0.5-8.0 keV) band images. The images are 
examined visually to locate other candidate sources, mainly close doubles and sources near 
the detection threshold. A total-band image reconstruction is generated using t he maximum 
likeli hood method to search for additional weak sources in the center region (jBroos et al. 



2010|). 



Refined sourc e positions and the final source list are established using the ACIS Extract 
software packag^l ( Broos et al. 2002 . 2010 ). Photons are extracted within polygonal contours 
of ~ 90% (~ 70% for crowded sources) encircled energy using position-dependent models 
of the point-spread function (PSF). Background extraction is performe d locally and take s 
into account spatial variations due to PSF wings of nearby point sourc es (jBroos et~aL 2010 ). 



Production of the final source list follows the simple procedure of iGetman et al] (12006 



20071 . l2008bl ) : the list of candidate sources is trimmed to omit sources with fewer than ~ 4 
estimated source background-subtracted counts (AC t /PSFfrac < 4; columns 7 and 11 in 
Tabled]). The final catalog has 225 sources. 

The ACIS Extract package is also used to construct source and background spectra, 
compute redistribution matrix files (RMFs) and auxiliary response files (ARFs), construct 
light curves and time-energy diagrams, perform a Kolmogorov-Smirnov (K-S) variability test, 
compute photometric properties, and perform automated spectral grouping and fitting. The 
following inferred source photometric properties are presented in Table [1] source coordinates, 
off-axis angle, net and background counts within the extraction region, source significance 
assuming Poisson statistics, effective exposure (correcting for telescope vignetting, satellite 
dithering, and bad CCD pixels), median energy after background subtraction, a variability 
indicator extracted from the one-sided K-S statistic, and occasional anomalies related to 
positions on chip gaps or field edges. Net counts has a range 4—1000 and typical values of 
~ 20. Median energy has a range of 1—5 keV and typical values around 2.4 keV. 



2 The A CIS Extract software package and User's Guide are available online 



at 



http : //www. astro .psu. edu/xray/docs/TARA/ ae_users_guide .html 
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3. Optical, Infrared, and Radio Counterparts 



3.1. Identifications 



Chandra source positions were compared with source positions from the 2-Micron All- 
Sky Survey (2MASS) NIR catalog (Figure [2]a,) . Chandra sources were considered to have 
stellar counterparts when the positional coincidences were better than 1" within ~ 3.5' 
of ACIS-I field center, and 2" in the outer regions of each ACIS-I field where Chandra's 
PSF deteriorates. Out of 225 X-ray sources, 190 have 2MASS matches. We have also 
acquired JHK images of the central part of the cluster from the United Kingdom Infra- Red 
Telescope (UKIRT), which have better spatial resolution and sensitivity than the 2MASS 
images (Figure [2b). Using UKIRT images 12 additional matches are identified, 10 of which 
are blended in 2 MASS. 



Two of the three OB stars identified in ISmith et al.l (119851 ) were matched to X-ray 
sources^. The third OB star was located on a chip gap and was not detected as an X-ray 
source. 



Rodriguez. Rodney. &: Reipurth I (120101 ) detect 20 Very Large Array (VLA) radio sources 
in the W40 complex near the center of the cluster, and 15 of them have X-ray counterparts. 
All sources with VLA cou nterparts also have infrared counterparts. The reported VLA source 
coordinates are rounded (IRodney fc Reipurthll2008l ). so Chandra-VLA positional offsets may 
be greater than 1" (#141). Individual radio and NIR counterparts for Chandra sources are 
listed in Table [2] along with their JHK S photometry. 

All X-ray sources with infrared counterparts are likely to be stars since YSOs have 
high apparent NIR fluxes compared to extragalactic contaminants. X-ray sources with NIR 
counterparts are classified as either W40 cluster stars or foreground stars (column 9 in Table 

Most of th e 23 X-ray detections without NIR counterparts are likely to be active galactic 
nucleii (AGN) (IGetman et al.ll2006l ). However, some of these are classified as W40 members. 
Sources #56, #140, #142, and #155 have strong variability, which is a characteristic of 
T-Tauri X-ray emission, and source #56 is located near the molecular core and is likely to 
be a highly embedded star. With their median energies > 3.5 keV, sources #56 and #155 



3 These OB sta rs are sources #99 and #141 in Table [TJ which correspond to OS 2a and OS la in 
Smith et al.l <|l985h . OS la is a blended source, where the components are distinguishable in the UKIRT 



image but not in the 2MASS image. The separation angle is 1.4" and the position angle is 14°. The position 
of the X-ray source corresponds to the position of the southern star. 
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are the most likely protostellar candidates, j udging from the evolu tionary class versus X-ray 
median energy relation shown in Figure 8 of iGetman et al.l ( 120071 ). 



3.2. IR Properties of W40 Stars 

Figure [3J presents NIR properties for most of the Chandra sources with NIR counter- 
parts based on the photometric data in Table [2j The nature of sources with known NIR 
photometry can be estimated from these diagrams. Our analysis uses models of intrinsic 
JHK S photometric properties for st ars with ages of 1 M yr a ssuming a distanc e of 600 pc 
(see SjEj). PMS evolutionary models of iBaraffe et al.l (119981 ) and lSiess et al.l (119971 ) are consid- 
ered for the mass ranges of 0.02 M < M < 1.4 M and 1.4 M < M < 7 M , respectively , 
and stellar photometric properties of stars with masses > 7 M are taken from ICoxl (120001 ) . 
Sources used in the analysis have errors on their NIR colors < 0.2 mag. 

The cluster stars are shifted to the upper right by Ay >5 mag on the color-color dia- 
gram (Figure [3^), so the foreground stars are easily identified due to the lack of significant 
absorption. Ten foreground stars are detected. 

The mass sensitivity limit of the X-ray detected PMS sample is ~ 0.2 M ; ~ 90% of 
stars have masses < 2 M ; and 7 stars have masses > 10 M© (Figure [3b). The 55 X-ray stars 
located >1.5 o to the right of the reddening vector originating at 0.2 M (Figure [3b.) are 
classified here as ii^-band excess stars^, leaving 109 X-ray stars as non i^-band excess stars. 
No stars with .K^-band excess are classified as protostellar objects; t he X-ray observation i s 
not sensitive to a protostellar population that might exist in W40 (IBontemps et al.ll2010l ). 
The disk fraction, taking into account the differences in X-ray sensitivity to stars with and 
without f^s-band excess, is estimated in 

Photometric mass and \^-band absorption derived from the color-magnitude diagram 
(Figure [3J) are tabulated in Table [2J Uncertainties in the photometric mass estimates are 
discussed in Appendix [XI Typical uncertainties are around 0.2 M for lower mass stars and 
0.8 M for higher mass stars. There is a degeneracy in JHK S photometry that leads to 
uncertain masses between 2 M and 10 M for the 1 Myr isochrone. We report the lower 
mass solution (2 M —4 M ) as masses in this range are more likely due to the shape of the 
IMF. 



4 Colors of massive stars identified as if s -band excess stars are located to the right of the reddening vector 
originating at the location of intrinsic colors of massive stars (this reddening vector is not shown in Figure 
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3.3. Membership 

Extragalactic sources, mostly AGNs at moderate redshifts, dominate X-ray source 
counts at high Galactic latitudes and some will be detected in fields near the Galactic plane 



despite the heavy obscuration (IBroos et al.l 120071 ). Even the low X-ray detection rate of 



main-sequence stars will lead to contaminants in the X-ray image due to the huge number of 
foreground field stars. The field stars detected by X-ray observations will tend to be younger 
members of the disk population, since stellar X-ra y emission declines rapidly after ~ 1 Gyr 



flFeieelson et al.l 12004 



Preibisch fc Feigelsonl 120051 ). However, background field stars are a 



less important source of contaminating X-rays for W40, due to the high obscuration from 
the cloud. 

Out of 225 X-ray sources detected by our observation 19 are classified as extragalactic 
sources ( §3.11) and 10 are classified as foreground stars ( §3.21) . The remaining 194 stars are 
likely W40 cluster members, although 30 do not have accurate IR colors for classification on 
the color-color diagram. Cluster membership is indicated in Table [2J 

The foreground and extragalactic contamination of the W40 sources may be compared 
with the analysis of contamination in a Chandra ACIS-I observation of the young-star cluster 
Cep B. The galactic latitude of Cep B and W40 are somewhat similar (~ 3.5° for W40 and 
~ 2.7° for Cep B), and the distances (600 pc for W40 and 700 pc for Cep B) and the Chandra 
exposures (38 ksec for W40 and 30 ksec for Cep B) are similar, but the galactic longitudes 
(28.8° for W40 and 110° for Cep B) are different. Confirmed by simulation s, the numbers 



of ex tragalactic and foreground objects in Cep B of 24 and 13, respectively (IGetman et al. 



20061 ). are similar to those of W40. 



YSC members are typically identified by IR excess, and K s -band excess is used to detect 
PMS stars within the ACIS-I field of view missed by Chandra. The 38 additional 2MASS 
cluster members with i^-band excess and the massive star OS 3a (no X s -band excess) are 
shown in Table [3J The i^-band excess is identified in a similar way as for X-ray sources, 
but a reddening vector originating from 0.09 M & instead of 0.2 M is used since the 2MASS 
point-source catalog is deeper than the Chandra observation and use of a reddening vector 
from 0.2 M would likely find foreground diskless stars. None of these 39 sources have a 
VLA counterpart. 



4. X-ray Spectral Analysis 



Low spectral resolution CCD spectra of young stars are typically modeled with one or 
two-temperature optically thin thermal plasma models subject to photoelectric absorption 
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to measure flux, temperature, and hydrogen column density. The spectral a nalysis of X-ra y 
stars is performed with the XSPEC spectral fitting package version 12.5 (lArnaudl Il996l ). 
The unbinned source and backgr ound spectra are fit with one-temperature plasma em ission 
models (apec, ISmith et al.ll200ll ) using the maximum likelihood method ( ICashl I1979T ) . We 
assume 0.3 times solar elemental abunda nces previously sugge s ted as typical for youn g stellar 
objects in other star formi ng regions (llmanishi et al.l l2001t iFeigelson et al.l 120021 ). Solar 
abundances are taken from 
the atomic cross sections of 



Morrison fc McCammonl (wabs, Il983h . 



Anders fc Gre vessel (119891) . X-r ay ab sorption is modeled using 



X-ray luminositi es and hydrogen col umn densities are also generated with XPHOT 
(XPHOT program in iGetman et al.l 120101 ) . XPHOT is a non-parametric method for esti- 
mation of apparent and intrinsic broad-band fluxes and absorbing X-ray column densities 
of X-ray sources. Measured X-ray luminosity and absorption values from both XSPEC and 
XPHOT methods are in good agreement. The advantage of the XPHOT method is that it 
is more accurate than forward-fitting methods for very faint sources and provides both sta- 
tistical and systematic (due to uncertainty in X-ray model) errors on derived intrinsic fluxes 
and absorptions. These errors on individual source X-ray luminosities are further used in our 
Monte-Carlo (MC) simulations to obtain errors on the X-ray luminosity functions of W40 
stars (SjHJ). 

The XSPEC and XPHOT spectral results are presented in Table HI The X-ray lumi- 
nosities included for objects classified as non cluster members may become useful in cases 
of source misclassifications. In addition, the XPHOT luminosities for OB stars may not 
be accurate because the XPHOT method assumes the T-Tauri X-ray emission mechanisms 
jGetman et alJboioh . 



5. Confirmatory Results: X-ray Luminosity, Stellar Mass, and Absorption 



We would like to improve our confidence in the multiwavelength study of the W40 
cluster by finding agreement between independently inferred X-ray and IR properties. Here 
the relationship between mass M and X-ray luminosity Lx and the relationship between 
visual extinction Ay and X-ray column density Nh are analyzed. 



The X-ray luminosity versus mass relationship was originally detected by IFeigelson et al. 

(119931 ) from ROSAT data. An empirical relationship of Lx oc M 17 , which extends over 3 



orders of magnitude in L x , can be seen in th e COUP observation of the Orion Nebula Cluster 
(IGetman et al.l 120051 ; iPreibisch et al.l 120051) and in the XMM-Newton Exten ded Survey of 
the Taurus Molecular Cloud (XEST) (iGudel et all 120071 : iTelleschi et al1l2007h . respectively. 
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This relationship has a poorly understood astrophysical cause, but may be an effect of the 
satur ation of the magnetic dynamo in the fully convective interior or on the surface of PMS 
stars (IPreibisch et al.ll2005l ). 



W40 cluster stars are plotted on top of the sample of lightly obscured PMS stars from 
COUP on the Lx-M diagram in Figure Hk- For the W40 T-Tauri stars the X-ray luminosities 
are the total-band (0.5-8 keV) luminosities generated by XPHOT (Q, and the masses are 
derived from the NIR color-magnitude diagram ( §3.21) . The W40 stars occupy the same 
parameter space on this graph as the COUP stars, down to the sensitivity limit of the 
Chandra observation of W40. 

IR absorption is sensitive to dust while bound-free absorption of X-ray photons mea- 
sures the integrated effects of interstellar material in any phase (partially ionized, neutral, 
molecular, solid). Contributions to the X-ray absorption come from ele ments including H e 
and inner shell electrons in C, N, O, Ne, Si, S, Mg, Ar, and Fe atoms (IWilms et al.l 120001 ). 
and the absorption is expressed as a hydrogen column density. The relation found be- 
tween Nh and Ay reveals the gas-to-dust ratio for a cluster. The relation between Nh and 
Ay in W40 (shown i n Figure HI) is in agreement with the gas-to-dust ratio of the galactic 
interstellar medium (IRyterl Il996l ). Both samples of stars with and without i^-band ex- 
cess follow the trend: log Njj/ log Ay = 21.30 ± 0.03 for stars with K s -band excess, and 
log Njj/ log Ay = 21.34 ± 0.02 for stars without fCj-band excess. Both the Lx — M and 
Nh — Ay agreements provide confidence in our X-ray and NIR analyses. 



XLF 



The statistical link between X-ray luminosity and stellar mass allows an association 
between the XLF and the IMF. A univer sal XLF has been hypothesized based on COUP 
and observations of IC 348 and NGC 1333 (IFeigelson et al.ll2005l ; iFeigelson fc Getman 1120051 ). 
he assumption o f a unive rsal XLF has been u sed to estim ate total po pulations of M17 
Broos et al.ll2007h . Cep B iGetman et al.lbood ). NGC 6357 JWangj I2007T ) . and NGC 2244 



(jWang et al 



20081 ). NGC 2244 and Cep B are complete beyond the COUP XLF's turnover 



at \ogL hc = 30.0 erg s" 1 (Figure E]), and while the NGC 2244 XLF closely matches the COUP 
XLF, there is a slight inconsistency between the Cep B XLF and the COUP XLF due to the 
excess of stars with log Lh c = 29.3 erg s -1 in Cep B. X-ray luminosity is known to decline 
during main-sequence evolution, likely causing the XLF of a cluster to change. However, The 
XLF does not strongly evolve during the PMS era for ages < 5 Myr (IPreibisch &: Feigelson 



2005 



Prisinzano et al.l 120081 ). The assumption of a uni versal XLF can also be used for 
YSC distance measurements (IFeigelson fc Getman II2005I ). For example the XLF analysis 
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for Serpens Main cluster by IWinston et al.l ( 120101 ) favors the distance of 360 pc over 2 60 pc, 
which is closer to the recent VLBA parallax measurement of 415 pc (IDzib et al.l 120101 ) . 



The W40 XLFs are simulated by varying reporte d XPHOT lum i nositie s of individual 
sources based on their errors, similar to the analysis of IGetman et al.l (120061). and are com- 
pared with the COUP XLFs of > 800 unobscured stars (IFeigelson et al.l 120051 ). Shifting the 
W40 XLF horizontally corresponds to changing the distance to W40, and the vertical shift 
of the COUP XLF model corresponds to scaling COUP to the W40 population. 



X-ray observations are less sensitive to stars with circumstellar disks (IGetman et al. 



20091 ). so we divide the stars into groups based on the presence of K s -baad excess; Figure 



[5] shows W40 XLFs of 50 stars with i^-band excess and 96 stars without i^ s -band excess 
in the hard and total bands. Both populations have similar XLF distributions above L/j C = 
10 30 - erg s" 1 (assuming D = 600 pc), but differ at lower X-ray luminosities. We also 
perform XLF analysis on the entire sample of 189 cluster members that have available X-ray 
luminosity information without dividing the sample based on K s -baiad excess (graph not 
shown) . 



6.1. Distance 

In Figure [5] the COUP XLF is compared to the W40 XLF with three different trial 
distances: 600 pc, the most commonly reported distance; 260 pc, the distance to the nearby 
on the sky Serpens South cluster; and 900 pc, the largest reported distance (JQ. Distances 
of 600 and 900 pc provide good fits, with a less complete 900 pc observed XLF. A distance 
of 500 pc provides a qualitatively worse fit, with a slight excess of stars above the COUP 
XLF at a log Lhc ~ 30 erg s _1 (graph not shown). With the assumption of a universal XLF, 
the W40 XLF analysis favors the distance of 600 — 900 pc; throughout the paper we adopt 
the most often reported distance of 600 pc. 

Further on we carry out the analysis assuming d = 600 pc. However, as mentioned in 
§1. H some studies favor the distance of 260 pc. Therefore, below we comment on possible 
modifications of the results if a distance of 260 pc is considered instead. 

If the distance to W40 is 260 pc, then the excess of stars with 28.5 < L hc < 30 erg s _1 
is inconsistent with a universal XLF. This may indicate the absence of the expected higher 
mass stars in the cluster (see §7] for a discussion of the IMF for 260 pc), or may be due to 
Lx evolution if the cluster is much older than 1 Myr. The excess wou ld be comparable to 



the excess of stars at Lh c = 29.3 erg s in Cep B ( IGetman et al.l 120061 ). 
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6.2. Completeness Limits, Total Population, Disk Fraction, and Age 



The X-ray completeness limits, i.e. where the W40 XLF diverges from the scaled COUP 
model, are log Lh c = 30.2 and logL te = 30.6 corresponding to 50% completeness at 1.4 M . 

The vertical shift between the COUP and the W40 XLFs for the total as well as the K s - 
band excess stratified star samples give the same estimate of the total W40 stellar population 
down to 0.1 M Q of 70% the COUP sample, i.e. 600 PMS stars (Table EJ. 

The same inferred vertical scaling of 35% to fit the COUP XLF model to both the 
W40 XLFs of stars with and without AVband excess gives a W40 AVband excess fraction 
of ~50%. The A s -band excess fraction found for the W40 cluster is likely to be a lower 
limit on the disk fraction, since JHK is not the most sensitive tracer of disks. For example, 
L-band photometry in addition to JHK gives a m ore robust estimate on the fraction of stars 
with circumstellar disks. Clusters such as IC 348 (IHaisch et al.lboOlbh ( age ~ 3 Myr), NGC 



2024 (IHaisch et al.ll2001af ) (age < 1 Myr), and the Orion Nebula Cluster (ONC; lLada et al. 
2000) (age 1-2 Myr) have the following reported JHK (JHKL) disk fractions: 20 - 25% 



(65%) in IC348, 60% (85%) in NGC 2024, and 55 - 90% (85%) in the ONC. The JHK S disk 
fraction of 50% thus suggests that the W40 cluster is young with an age of likely < 1 Myr. 
W40's likely astr o physi cal disk fraction of ~80% gives an age of 0.5 Myr using the regression 
line in iMamajek! (120091 ) for age versus disk fraction of 22 star clusters. 



However, there is a somewhat controversial situation, since 6 out of the 8 high-mass 
star candidates (M > 10 M Q ) have A s -band excess, suggesting a A s -band excess fraction 
of 75^^ 9 %0. Theoretical disk depletion timescales for early-type stars, due mostly to photo- 
evaporation with some contribution from material removal by winds, are rough ly 0.2 Myr 
and 0.7 Myr for 07V and B1V stars, respectively jHollenbach et al.lll994[ l200oh . The ob- 
served disk fractions and the theoretical depletion timescales might suggest that most of W40 
X-ray star high mass candidates formed somewh at later than low mass stars ; for example 
the latter phenomenon is reported for stars in W3 
Povich fc Whitnev! (hoioh . and IRAS 19343+2026 



Feieelson & Townslevi (120081 ) . M17 SWex 



Oiha et all f]2010h . 



If W40 were at 900 pc the XLF would indicate completeness limits at log L^ c = 30.4 erg s" 1 
and logL te = 30.8 erg s _1 and a total population of at least 1000 stars down to 0.1 M Q , 
but the disk fraction would not change. However, if the distance is 260 pc then these values 
cannot be derived from the XLF since the shape of the XLF would differ from the shape of 
the COUP XLF. 



5 Errors on disk fra ctions have been estimated using binomial distribution statistics as described by 
Burgasser et al.l (|2003l ). 
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IMF 



To verify the results of the XLF analysis, the IMF of the cluster is inspected (assuming 
a distance to W40 of 600 pc). For consistent treatment, in the IMF analysis the stellar 
masses used for both W40 and COUP are the masses derived by the method described in 
§3.21 The W40 X-ray and IR selected s amples are com bined and compared to the sample 
of > 800 unobscured COUP stars. The Ichabrier! J200J IMF is scaled to the COUP IMF, 
the IMF of sources without K s -band excess, and the IMF of sources with i\~ s -band excess 
(Figure Ekb). Error bars on the W40 IMF are generated through simulation^. 

W40 cluster members with no i\~ s -band are complete to 0.5 M and stars with i^-band 
excess are complete to 1.5 M Q . This agree s with the completeness result from the XLF 
analysis. From the scaling of the IChabrierl (120031 ) IMF, the population of stars without 
K s -ba.ri.d excess (with i^-band excess) is 32% (40%) the size of the COUP sample, implying 
a JHK S disk fraction of 55% and a total population of ~600 PMS stars down to 0.1 M Q 
(Table E}. These results are consistent with those of the XLF analysis. 

If W40 is at a distance of 260 pc, the IMF for both stars with and without i^ s -band 
excess would be complete down to 0.1 M , so the total cluster population would be ~200 
stars (Figure Ebd). The disk fraction would be ~50%. However, at this distance the IMF for 
stars without i^-band excess has an excess of stars between 0.1 and 0.3 M , so the IMF does 



not resemble either the COUP IMF or the IChabrierl (120031 ) IMF. This result also indicates 
that 260 pc is too near a distance for W40. 



8. Massive Star Candidates 



W40 is known to contain massive st ars, including three OB stars that act as ionizing 
sources for the cluster (ISmith et al.lll985l ). Eight cluster members have photometric masses 
> 10 M (Table El), and the masses of identified ionizing sources OS la, OS 2a, and OS 3a are 
consistent with OB stars. OS la is a visual double and only binary IR photometric properties 



6 The simulation takes into consideration effects of statistical errors, systematic errors, and mass degen- 
eracy for some stars. The selection of stars in each mass bin is generated 1000 times with individual mass 
values calculated from JHK S photometry where magnitudes are drawn from a Gaussian distribution with 
a standard deviation equal to error reported by 2MASS. For stars with M < 5 M©, masses are drawn from 
a Gaussian distribution with a mean equal to the mass derived above and a standard deviation of 30% 
to account for systematic error found in Appendix Rl For star s with multiple mass solutions, a solution 
is randomly chosen with a weight based on the (|Chabrierl 120031 ) IMF. However, the effect on the IMF is 
minimized since a single mass bin spans the entire region of degeneracy. 
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are known. OS 2a has a mass of 30 M and OS 3a (which is not observed in the X-ray due 



to location on the chip) has a mass of 15 M . Out of the 8 hig 
have .fTs-band excess, including OS la and OS 2a, but not OS 3a. 



i mass star candidates, 6 



Smith et al. (1985) detect 



circumstellar material around OS la, OS 2a, and OS 3a, and ( IVallee fc MacLeodlll994l ) find 
evidence for substantial circumstellar material around OS la and OS 2a. 

However, if a distance of 260 pc is considered only two stars, OS la and OS 2a, have 
JHK S masses > 10 M . 

X-ray properties inferred from one-temperature spectral fita3 are also given in Table |6j 
We caution that due to the high absorption of the W40 cluster, it is difficult to infer the 
existence of very soft temperature components. In addition lack of low energy photons may 
bias temperatures inferred by XSPEC to higher values. 

Soft (< 1 keV) and nearly constant X-ray emission from O and early-type B stars is 
often explained by a wind-shock m odel, where shocks form in a freely expandin g wind sub- 
ject to the Line-Driven Instability ( jLucy fc White! Il980l ; lOwocki fc Cohenlll999l ). Hard and 
often variable X-ray emission from massive stars is propose d to be produced by eit her the 



wind-wind collision mechanism in a massive binary system ( Zhekov fc Skinner 



1997 



2000h or the 



Gagne et al. 



Magn etically Channeled Wind Shock Mechanism ( IBabel fc Montmerle 
20051 ). None of the high mass star candidates is X-ray variable. Considering large uncertain- 
ties of ~ 50% on inferred plasma temperatures, spectra of two sources, #122 and #145, are 
significantly inconsistent with the purely soft (< 1 keV) plasmas of the classical Line-Driven 
Instability model. Highly ionized spectral lines of Mg, Si, and S with formation temperatures 
above 2 keV are also seen in the spectrum of source #122, suggesting the presence of the hot 
plasma component. An alternative to the wind- wind collision or Magnetically Channeled 
Wind Shock Mechanism explanations is that the hard X-ray emission is not coming from 
the OB stars themselves, but from unresolved T-Tauri companions. 

The level of X-ray emission of O and early-type B stars is often com parable to or 
higher than the most X-ray luminous T-Tauri sta rs (e.g. Stelzer et al ^ 
emiss ion roughly follows the trend L x = 10~ 7 L Bo i ( IChlebowski et al.lll989t iBerghoefer et al. 



19971 ). Figure [7] shows the L tc — L bol relation for the W40 T-Tauri stars and OB candidates. 
Bolometric luminosities are derived from the NIR color-magnitude diagram ( §3.2p .The mean 
log(L te /LBoi) for W40 OB candidates is -7.3 and the distribution has a standard deviation 



7 Results are given for one -tempe r ature thermal fits (wabs * apec). To account for mild pile- up in source 
#122 the model pileup from IPavisI ( 2001 ) is applied. Reduced x 2 is used to fit models to all the so urces, 
excep t for sources #46 and #144 which have a small number of counts and are fit using the C-statistic (jCash 
1979) on the unbinned spectrum. 
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of 0.7. This i s similar to the re lation found for other clusters obs erved by Chandra, for 
example M17 faroos et al.l hpOTl ) and NGC 2244 f lWang et all l2008h . However, the X-ray 
emission level of some W40 high-mass candidates may be explained by an X-ray active 
T-Tauri companion unresolved in the NIR. 



9. Spatial Structure 
9.1. General Cluster Properties 



For the analysis of W40's structure, we choose to define the center of the cluster as 
the median position of all C/iandra-detected members, which has coordinates R.A.=18 h 31 m 
25fl4, decl.=-2° 05' 35"7 (J2000). An alternative strategy is to use the OB stars as the 
center of the cluster, but in W40 the median position of the high mass stars is offset from 
the median positions of other stellar populations as shown in Table [7J so we judge that 
the median position of all the stars is a better estimate. For discussion of projected radial 
distances in the cluster d = (-D/600 pc) accounts for uncertainty in the distance D to W40. 

Fifty percent of the C/iandra-detected cluster members are within 0.50c? pc of the cluster 
center and 75% are within 0.90d pc. The edge of the ACIS-I array is 7.5' (1.3d pc) from the 
center of the cluster, so our analysis of structure is restricted to the region within this radius. 
The most distant detected star is 1.8d pc from the center. The core radius is ~0.15— 0.2d pc 
In the core there is a small clump of 8 stars, including OS la, within a radius of 10". 



A smoothed stellar surface density plot of the cluster is shown in Figure [HJ 

The structure of YSCs may be affected by both the locations where stars form and by 
cluster dynami cal evolution. The v elocity dispersion of stars in W40 is likely to be similar 



to other YSCs; iFuresz et al.l (120081 ) find a radial velocity distribution of a = 3.1 km s 1 in 



the ONC, and simulations of small and large clusters find one- dime nsional stellar velocity 
dispersions of 1.2 km s _1 and 4 km s _1 (IBate et al.ll2003l ; lBatdl2009l ). Assuming the stellar 
velocity dispersion in W40 of ~2-4 km s _1 and a cluster size of 2d pc, the cluster crossing 
time is t cross ~ 0.5 — Id Myr, on the order of the age of the cluster. The two-body relaxation 
time for a cluster is 



t 



N 



relax 



81niV 



X t r 



where N is the number of stars ( iBinney fc Tremaind 119871 ) . In W40, N~600 so the two- 
body relaxation time is ~ 5 — 10 Myr , longer than the a ge of the cluster. However, violent 
relaxation may proceed more quickly ( iLynden-Belll 119671 ) . 



If the distance of 260 pc is considered instead of 600 pc, the relaxation time would 
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decrease to ~ 1 — 2 Myr. 



9.2. Surface Density Profiles 

9.2.1. Radial Profile 



Figure M shows the radial surface density profile of W40. All cluster members detected 
by Cha ndra are included to improve counting statistics. The radial profile of stars from 
COUP (IFeigelson et al.ll2005l ) is provided for comparison, and a scaled version of the COUP 
radial profile closely resembles the radial profile of W40. Various functional profile models 
are evaluated using the cumulative distribution function (CDF) of the radial distances of 
cluster members and the K-S statistic. No single power law provides a good fit. However, 
the inner 0.2c? pc part of the profile is flatter and can be fit with a power law index of -0.2, 
and the outer part of the cluster is fit with a power law index of -1.6. The 90% confidence 
interval on the indices is ±0.3. YSC surface-density profiles are often fit by shallower power 
laws that typicall y have an exponent of -1. Examples in clude IC 348 (a oc r _1 ) and NGC 
2282 (ffocr 11 ) (ILada & Ladalll995l : EoTner et al.lll997h . 



Assuming spherical symmetry, a surface density distribution with a power law a cx r a 
and a < will correspond to a three-dimensional density distribution with a power law 



p oc r 



a-l 



so the three-dimensional power law index for the inner part of the cluster will 



be approximately -1.2 and the power law index for the outer part of the cluster will be -2.6 
(similar to the three-dimensional power-law index derived from the Q parameter analysis 
described in 391T2]) . 

We define the core of W40 as the region within 0.2c? pc, since this region has a noticeably 
less steep profile than the outside of the cluster. Within the core there are 48 Chandra- 
detected stars, the three ionizing sources OSla, OS2a, OS3a, and 3 other high mass star 
candidates. This core radius is similar to the core radius of the ONC which was found to be 



0.15 pc ( [Hillenbrand fc Hartmann! 119981 ). 



The ONC stellar density pro file is well fit by the King profile ( IHillenbrand & Hartmann 



19981 ) for an isothermal sphere ( IKing I Il962f ) that is applicable to gravitationally relaxed 
clusters and is described by the equation 



a r 



OAT] 



1 + (r/r c y 



(2) 



Nevertheless, IHillenbrand fc Hartmann! (119981 ) , IFeigelson et all (120051 ) , and lFuresz et all ( 120081 ) 
argue that the Orion cluster is not truely dynamically relaxed. The W40 cluster can also 
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be fit by an isothermal sphere with a core radius of 0.15d pc (K-S probability of 92%), but 
W40 might also not be dynamically relaxed due to its young age ( §9.11) . 



9.2.2. Profile Fitting and Detection of Subclustering with the Q Statistic 



A statistic Q indicates the existence of subclustering and measures t he central concen- 
tration of clusters without subclustering (ICartwright fc Whitworthll2004l ). This parameter 
combines a measure of the separation between stars and a measure of the edge lengths of the 
Euclidean minimal spanning tree (EMST). The EMST is the unique graph with minimum 
total edge length that connects all vertices. Figure [10 shows the EMST for all W40 cluster 
members detected by Chandra, where the stars are the vertices of the graph. The param- 
eter s is the mean separa t ion b etween stars, and m is the mean edge length on the graph. 
Cartwright fc Whitworthl (120041 ) define 



s/r B 



(3) 



and 



m 



m(N n 



N„ 



(4) 



where r max is the maximum radius, which is 7.5' for W40 (stars outside of this radius are 



excluded), and N n 
by 



is the number of sources within the radius. The parameter Q is defined 



s 



(5) 



Figure 5 in ICartwright fc Whitworthl ( 120041 ) gives fractal dimension or central concentration 
as a function of Q generated from simulated clusters. A YSC with subclustering will have 
a value of Q < 0.8. For W40 Q = 1.47, indicating that W40 does not have subclustering 
(within 7.5') and that the cluster has a three-dimensional density profile p oc r a , where 
a = —2.7. In addition, a visual examination of the EMST in Fig [10] reveals no unusually 
long segments, which suggests that no subclustering is present. The analysis also agrees 
closely with the results of the power law fits to the outer part of the cluster ( §9.2.ip . For 
comparison, the shallower three-dimensio nal radial profile exponents of -1-2, -2.2, -1.7 ar e 
found for p Ophiuchus, IC348, and a Ori JCartwright fc Whitworthl I2004J : ICaballeroll2008l ). 
respectively. 
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9.2.3. Cluster Elongation and Asymmetry 



There is no evidence of that the W40 cluster of all CTiemdra-detected PMS stars is 
elongated; no correlation is found between right ascension and declination of cluster members 
(Person correlation coefficient = -0.1), and the standard deviations of right ascensions and 
declinations are similar (173" and 180", respectively). 

To further investigate whether the exterior of W40 (r > 0.2c? pc) has radial symmetry, 
the distribution of azimuthal angles of stars is compared to an isotropic distribution, and 
the median distances of sources from the center is examined as a function of azimuthal angle 
(Figure ITT|) . For the former test the CDF of the angles wa s found to be consistent with a 
uniform distribution using Kuiper's statistic (IKuiper Ill962[ ). providing additional evidence 
that subclustering does not exist in W40. For the latter test, sources are binned by azimuthal 
angle with bin sizes of 60°, and median distances are determined with error on the median 
calculated using 1000 bootstrap values. The median distances are consistent with a single 
value (reduced y 2 is 0.4 when fit with a constant). However, with bin sizes of 180° the 
reduced x 2 is 3.4 due to larger radial distances in the south half of the cluster, indicating 
that these stars are less centrally concentrated. 



Most YSCs, including both the ONC ( Hillenbrand fc Hartmannl Il998l ) and o Orionis 



Hillenbrand & Hartmann 


1998) 


g molecular clouds ( 


Allen et al. 



20071 ). If W40 is not dynamically relaxed, the near radial symmetry could be a result 

pi 

■ ( 



of spherically symmetric star for mation or a project ion effect if we are looking down the 
asymmetric axis of the cluster. iFuresz et al.l ( 20081 ) caution about the interpretation of 



radial profiles that average over all angles, since asymmetry may cause steeper power laws, 
but this is not a great concern for W40 because of the near radial symmetry. 

Finally, if considering only sources with i^-band excess the azimuthal distribution is 
not uniform (0.5% K-S probability). We further study this in §9.41 



9.3. Mass Segregation 

The dependence of cluster structure on stellar mass in YSCs provides clues about clus- 
ter formation and dynamical evolution. It is observed that massive stars are preferentially 
detected near the centers of open clusters and YSCs. Mass segregation is an expected 
result of dynamical relaxation, but must be explained by different astrophysics for its ap- 
pearance in clusters that are not dynamically relaxed. Mass segregation is seen in clus- 
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ters that are not dynamically relaxed, such as NGC 2071, the Trapezium , and NGC 2074 
( ILada et al.l Il991t iHillenbrand fc Hartmannl Il998l ; iBonnell fc Daviesl Il998l ) . This mass seg- 
regation would be caused i f higher mass stars tend to be formed nearer the center of clusters 
rtProszkow fc Adamsl 120091 ) or it could be an effect of dy namics during the f o rmati on of a 
not-yet dynamically relaxed star cluster as described by iMaschberger et al.l (120101 ). This 
phenomenon is also seen in W40. 

We stratify W40 cluster members by 3 mass ranges at 1.5 M and 10 M and by K s - 
band excess (we will refer to the mass ranges as low mass, intermediate mass, and high mass). 
Both the intermediate-mass and high-mass populations should be approximately complete 
(based on the results of the XLF/IMF analysis). The spatial positions of stars in these strata 
are plotted in Figure [T2] and the centers and median radii are given in Table UJ Low-mass 
stars are clearly less centrally concentrated than medium-mass stars and high-mass stars. 
The high-mass stars are the tightest group, except for star $46 which is several arcminutes 
away from the rest of the group. We do not have proper motions to evaluate whether this 
star is a runaway. 

The CDFs of radial distances of low-mass stars, intermediate-mass stars, and high-mass 
stars show these mass segregation trends, and the K-S probab ilities that the thre e mass 
distributions come from the same population are negligible^. Moeckel fc Bonnell ( 2009 ) 
describe a method that uses the area of the convex hull for n stars to search for mass 
segregation. The convex hull is defined as the minimal convex polygon which covers a set of 
points, and the area may be used to determine stellar density. For W40, the probabilities 
that n randomly selected low-mass stars could have a convex hull with an area less than the 
area of the convex hull of the n high or intermediate-mass stars is 7 x 10~ 4 and 5 x 10~ 2 , 
respectively. For the 25 stars with M=1.5±0.1 M Q the probability is 2 x 10~ 2 . 

This result is unusua l since mass stratific ation typically is seen starting at higher masses 
than 1.5 M , for example lAUison et al.l ( 120091 ) find no mass segregation in the ONC at masses 
less than 5 M^. However, mass segr egation beginning at 1 M Q in 1 Myr old simulated clusters 
is found by iMoeckel &: Bate! (120101 ) in an N-body simulation. 



The distance to W40 does not affect the existence of mass segregation, since radial 



8 Each of the samples (low mass, intermediate mass, and high mass) is generated 1000 times using the 
same method described in Sj7]to take into account statistical and systematic errors on mass. The K-S test 
of the radial-distance CDFs is repeated for each of the 1000 simulations, and we find that more than 90% 
of the significance levels are negligible for comparison of the low mass sample to the intermediate and high 
mass samples. The intermediate and high mass samples have marginally different CDFs in the simulated 
data. 
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distances scale similarly and the masses derived for 260 pc are an increasing function of the 
masses derived for 600 pc. However, mass segregation would starts at 0.3 M for 260 pc 
rather than 1.5 M . 



9.4. Distribution of i^-Excess Stars 

In Figure [T2a the low mass stars with K s -band excess are concentrated on the east 
side of the cluster, while the stars without i^-band excess are evenly spread over the entire 
cluster. This effect is not as strong among the intermediate mass stars, many of which reside 
in the core from both populations. In Table [7] the median position for low-mass stars with 
i^-band excess is 0.27c? pc to the east of the cluster center, which is farther away from the 
center than for any of the other populations listed. 

Figure [12] shows the spatial distribution of K s -excess stars and the molecular core. 
Almost no stars with .fC,-band excess are located west of the cluster core, and 20% of the 
stars without i^ s -band excess are located west of the westernmost i^ s -excess star. The K-S 
probability that the spatial distributions are the same is P = 2 x 10~ 6 . 

To the west of right ascension 18 h 31 m 15 s (277.81°), 23 Chandra-detected stars have no 
i^-band excess and Chandra-detected stars have i^-band excess, s o the la uppe r limit on 



the observed JHK S disk fraction in this region is 5% (equation 21 in I Gehreld 1 1 98 6l ) . East of 
18 h 31 m 24 s (277.85°), 45 Chandra- detected stars have no K s -band excess and 35 Chandra- 
detected stars have -K^-band excess, so the observed JHK S disk fraction in this region is 
~ 44^7%. Since only 20% of T-Tauri stars with ii^-band excess are detected, compared to 
45% of T-Tauri stars without ii^-band excess, the astrophysical JHK S disk fraction in these 
two regions is < 10% and 63%, for west and east regions respectively. 

The extremely different disk fractions must have been caused recently, since this phe- 
nomenon will be erased due to mixing on the order of the cluster crossing time, which is 
< 1 Myr. This provides further evidence that the age of the population of stars with i^-band 
excess is <1 Myr. 

Different -fT s -band excess fractions may be attributed to environments hostile to disks 
or to different stellar ages. In the second case, our assumption of a coeval cluster is violated. 
This may be due to two different epochs of star formation or progressive star formation 
moving from west to east. Alternatively, the spatial distribution could be explained if the 
youngest stars moved to the east side of the cluster due to a similar initial velocity. Sub- 
clustering is also a possibility, although §9.21 gives evidence that the stars are from only one 
cluster. 
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A depleted disk fraction due to a harsh environment in the western part of the cluster 
seems unlikely, since the i^-band excess disk fraction is high in the cluster core where the 
UV radiation and stellar winds would be strongest, and even some of the most massive 
stars have i^-band excess. The disky stars in the core could have moved there recently, so 
they would not have been exposed to the environment for long enough to lose their disks. 
However, many of the disky stars in the core are high or intermediate mass stars, which are 
concentrated toward the center of the cluster ( §9.3j) . 



9.5. Dust Lane 



in 



The MSX data show th at W40 is a bipolar nebula with an hourglass shape (Figure 6 



Rodney fc Reipurthl 120081 ). This type of structure is fairly c ommon for H II regions; 



for example W40's shape close ly resembles Sharpless 106 (S106) (lOasa et al.ll2006l ). IRAS 
19343+2026 Joiha et alihoioh . and the Lagoon Nebula foothill et al.ll2008h . In W40, the 
waist of the hourglass is the location of a dust lane which bifurcates the H II region and 
shows up in the optical image (Figure [T4"b). Larger column densities of dust and gas are 
also measured for stars in the dark lane from the JHK S data and the X-ray data. The 
arrangement of the J — H >2.2 sources (corresponding roughly to Ay > 10 — 13 mag) is 
shown in Figure [T4h in comparison to the rest of the cluster and to the molecular core. We 
draw a region which encloses most of the J — H >2.2 sources, but few of the less absorbed 
sources, and plot this on a map of the smoothed extinction.] (shown in Figure ITS]) . The 
mean XPHOT hydrogen column density is (3.2 ± 0.4) x 10 22 cm -2 in the dark filament and 
(1.8 ± 0.1) x 10 22 cm -2 outside this region, considering only sources without i^-band excess. 
These correspond to visual extinctions of 15 mag and 9 mag, respectively. The filament 
appears to have a normal gas-to-dust ratio: log Nh/ log Ay = 21.33 ± 0.03 for stars in the 
filament, and log Nh/ log Ay = 21.34 ± 0.02 for stars outside the filament. 

The dust lane passes through the eastern half of the cluster core, and lies several arcmin- 
utes to the east of the CO core. The filament has dimensions 3' by 12', and an approximate 
angular area of 36 arcmin 2 or 1.1 xlO 37 cm 2 at a distance of 600 pc. The filament's mass is 
estimated by multiplying its angular area by the difference in mean column density inside 
and outside the lane giving a total mass of ^200 M . This is comparable to the mass of 
the molecular core found by IZhu et al.l (120061 ) . If we assume the filament's extent along the 
line of sight is similar to its width of 3 arcmin, then the hydrogen density of the filament is 



9 The map of Ay in the cluster is generated by smoothing the Ay of stars with a Gaussian kernel with 1' 
standard deviation. 
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10 4 atoms cm 3 . 



If the distance of 260 pc were considered instead, the masses of both the molecular core 
and dust lane would be scaled down by a factor of ~5, but the density would be the same. 

The presence of a large number of sources with K s excess near the dust lane, suggests 
that the lane may be a region of active star formation. This is f urther supported by the high 
concentration of Herschel protostar candidates from Figure 3 of iBontemps et al.l (120101 ) that 
lie in this region. In contrast, few deeply embedded X-ray sources are found in the CO core, 
suggesting that the molecular core is not the location of present star formation. 

The geometry of the dust la ne and the molecular core sug gests that they form a ring, 
similar to the rings described by iBeaumont &: Williams! (120101 ) . However, YSOs projected 
on the molecular core do not have very high Ay, and CO is not detected strongly in the 
dust lane. The former problem may be explained if the ring is oriented from our point of 
view so that the dust lane passes in front of the cluster and the molecular core is behind the 
cluster. The molecular core probably has a normal gas-to-dust ratio and we would see both 
CO emission and a visual dust lane if we were observing the cluster along a line of sight that 
passes through the cloud. 

The lack of CO detection in the dust lane might mean that the molecular material 
has been destroyed in this region, or it might be an effect of the lack of sensitivity of 
the CO observations. CO may become froz e n out onto dust grains in a molecular cloud 



(Ivan Dishoeck fc Blakd 119981 ). iMitchell et al.l (120011 ) find evidence of this effect in two parts 
of the Orion B molecular cloud with bright dust continuum but almost absent CO emission. 
Alternatively, CO may be destroyed by UV radiation or expanding HII regions produced by 
OB stars. 



10. Interesting X-ray Sources 

10.1. Stacked Point-Source Spectra 

Out of 225 X-ray sources, only 24 sources have more than 100 net counts, and almost 
a third of the sources have fewer than 10 net counts. Composite spectra are useful for 
gaining insight into the sources which lack sufficient counts to be fit parametrically (e.g. 



Feigelson et al.ll2005l ). They are created for the W40 point sources by merging individual 
source spectra, backgrounds, RMFs, and ARFs using the ACIS Extract software. Separate 
composite spectra are generated for cluster stars with more than 100 counts, cluster stars 
with fewer than 100 counts, and the extragalactic sources; these spectra are plotted in Figure 
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H6b. 



The spectra are fit in XSPEC using the thermal-plasma model apec and the inter- 
stellar photoelectric-absorption model wabs. A two-temperature model with fixed parame- 
ters Z = 0.3 Z ro (see §3} and kTi=0.8 keV (typical cool component in T-Tauri stars from 



Preibisch et al.l ( 120051 )) best fits the composite spectra of >100-count (<100-count) sources 
with N H = 1.7 x 10 22 cm" 2 and kT 2 =5.4 keV (N H = 2.2 x 10 22 cm" 2 and kT 2 =5.2 keV). 
For >100 count (<100 count) sources, model fitting with metal abundance Z as a free pa- 
rameter gives a value of Z =0.43+g;Qg (Z =0.44^;^), indicating that Z=0.3 Z is a good 
approximation for the cluster. The stacked spectrum of extragalactic sources is well fit by 
an absorbed po wer law with a photon index of 1, within the typical range of 1 — 2 for active 
galactic nuclei ([Alexander et al.l 120031 ). 



For the composite spectrum of >100-count (<100-count) sources XSPEC calculates ob- 
served (not corrected for absorption) hard-band fluxes of (2.9±0.1) x 10 -12 erg s _1 cm -2 
((1.7±0.1)xl0~ 12 erg s -1 cm -2 ) and soft-band fluxes of (3.5±0.2) x 10~ 13 erg s _1 cm -2 
((2.6±0.1)xl0~ 13 erg s _1 cm -2 ), implying observed hard-to-soft flux ratios of 8. 3±0. 6 (6.5±0.5). 



10.2. The Diffuse X-ray Emission 



In a number of star forming regions the soft (kT < 1 keV) diffuse thermal emission has 
been detected from a gas shocked in w ind- wind or wind-surrounding medium interactions 
( iTownsley et al.ll2003t iGiidel et al.ll2008l). The uncommon hard diffuse synchrotron emission 
has also been reported (e.g. IWolk et al.l 120021 ). In some cases the observed X-ray diffuse 



emiss ion can be completely explained by the unresolved stellar population (e.g. lGetman et al. 



20061 ). 



The Chandra observation W40 has detected diffuse X-ray emission near the core of 
the cluster that is unassociated with the resolved point sources and can be seen in the 
smoothe d image in Figur e |5J To isolate the diffuse emission, a "Swiss-cheese" mask is 
created ( jBroos et al.ll2010l ). Despite the Swiss-cheese masking, there is some light remaining 
from the wings of very bright X-ray sources concentrated in the core. To avoid any of this 
contamination a generously-sized polygon surrounding the central stars is excluded (Figure 
[T6"a). Three background regions are used; backgrounds 1 and 2 are on the the Chandra ACIS- 
13 chip, while background 3 is on the Chandra ACIS-I2 chip. Ba ckg round 2 is adjacent to 
the extraction region, while backgrounds 1 and 3 are further awa 



ackgr 

S 



'In the smoothed image with all point sources removed, there is a hint of extra diffuse emission to the 
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A moderate- resolution spectrum is extracted with ACIS Extract from the source region 
and the three background regions. Fits with different backgrounds give similar results; 
reported results are given for background #3. A degeneracy in spectral model fit exists due 
to the high absorption which makes it difficult to constrain the soft component (kT < 1 keV) 
due to the small number of counts. A one-temperature model produces a good fit (reduced 
X 2 =0.94) with N H = (1.2 ± 0.2) x 10 22 cm" 2 , kT= 2.9 ± 1.0, and Z = 0.3 Z Q . A variety of 
two-temperature models work as well, including Njj = 2.2 x 10 22 cm" 2 , kT\ = 0.66, kT 2 > 15, 
and Z = 0.3 Z and N H = 3.3 x 10 22 cm" 2 , kTi = 0.27, kT 2 = 4.3, and Z = 0.3 Z , both 
with a reduced x 2 °f ~ 0.8. 

The observed (uncorrected for absorption) hard-band and soft-band fluxes of the diffuse 
emission are (2.0 ± 0.2) x 10~ 13 erg s _1 and (4.8 ± 0.1) x 10 -14 erg s" 1 , respectively, so the 
hard-to-soft flux ratio is 4.2 ± 0.4. This is smaller than that of the stacked stellar X-ray 
sources at a statistically significant level (P < 0.001), which demonstrates that the spectra 
have a different shape. This can be seen visually in Figure [16b, where the bump at low 
energies is relatively larger in the diffuse spectrum than in the stacked stellar spectra. This 
result is an indication that the diffuse spectrum may not be entirely stellar. However, the 
unresolved stellar population, being composed of lower mass source s, is also expected to have 



spectra that are softer than the resolved population (Figure 11 in lPreibisch et all 120051 ) 



In the hard band only X-ray emission from stellar sources should contribute to the diffuse 
emission. The measured hard-band luminosity of the diffuse emission is compared with the 
expected hard-band luminosity of undetected stellar sources within the extraction region. 
11% of the C/iandret-detected stars lie within the extraction region, so we assume a similar 
percentage of the unresolved stars lie in this region. The sum of the hard-band luminosity 
of all the missing sources from the XLF below the completeness limits is 6.8xl0 31 erg s _1 , 
and thus the hard-band emission expected within the extraction region is 7.5 xlO 30 erg s _1 . 
The hard-band luminosity of the diffuse emission is calculated using the one-temperature and 
two-temperature models discussed above, and uncertainty is estimated by varying parameters 
within their 1 a errors. The result is Lh c = (8.2 ± 3.4) x 10 30 erg s _1 , which is consistent 
with that of the unresolved stellar component. 

The estimated total-band luminosity of undetected T-Tauri stars in the extraction region 
is 2.0 x 10 31 erg s" 1 , which is calculated in the same way as the estimated hard-band lumi- 
nosity. The total-band luminosity of the diffuse emission is L tc = (1.9 ± 0.7) x 10 31 erg s _1 



south- west of the main diffuse emission. These excess counts lie near sources #48 and #51, within the 
molecular core. A possible origin could be a group of highly embedded PMS stars. However, spectral 
analysis suggests that the source spectrum is not significantly different from the local background. 



-26 - 



for the one-temperature model fit, consistent with that of the undetected T-Tauri stars. 
However, it may be more than five times greater for the two-temperature model fits. Al- 
though, the total-band diffuse emission may be consistent with only the emission expected 
from unresolved stars, it is unclear if a real nonstellar X-ray component is present due to 
the degeneracy in spectral model fits. 

If the distance to W40 is 260 pc, then diffuse emission luminosities decrease by a factor 
of 5. However, the IMF analysis indicates that the observation would be complete down to 
0.1 M Q so the X-ray emission from unresolved stars would be an order of magnitude lower 
in both the hard and the total bands. 



10.3. Super-Hot Flare in Chandra Star #138 



Long and powerful flares from YSOs are sometimes detected in Chandra observations of 
YSCs; for example, in the COUP proj ect ~200 flares are de tected indicating a flare rate of 
approximately one per star per week (IGetman et al.ll2008al ). During the 38 ks observation 
of W40, the decay phase of a powerful flare is detected in the light curve of source #138, 
an intermediate mass (4 M Q ) YSO near the center of the cluster. This flare has a peak 
luminosity >40 times its non-flaring luminosity and an e-folding decay timescale of 7.6 ks. 
The flare is strongly affected by pile-up on the ACIS-I chip, so we do not present flare cooling 
analysis. 

At least two alternative explanations for powerful flare activity from a 4 M YSO 
exist. The object can be an F-type T-Tauri star with a partially developed radiative core. 
Big super-hot flares are typically seen from two types of T-Tauri stars: (1) low- mass stars 
(M < 2 Mfp) with a ctive disks, (2) and more massive (M > 2 M Q ) disk-free stars (§5 in 
Getman et al.ll2008a ). It is possible that the formation of a radiative core in the latter case 
is associated with the change in the magnetic field dynamo, for example from turbulent 
to alpha-omega, which in turn strengthens their magnetic fields and allows production of 
super-hot flares. 

Alternatively, the object can be a fully radiative A or late B-type star. Although X-ray 
activity in late B and early A stars is generally wea k, X-ray flares asso ciated with these stars 
have b een seen in several cas es such as HD 161084 (IMiura et al.ll2008l ). HD 261902, and HD 
47777 (jYanagida et al.l 120071 ). One of the possibilities for the X-ray emission is the presence 
of an X-ray active T-Tauri companion to an IR-bright, A/B-type star. 
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11. CONCLUSIONS 

We present a 38 ksec Chandra ACIS observation of the W40 complex. Although W40 
is one of the nearest massive star-forming regions, it is poorly studied due to its high obscu- 
ration. Our main conclusions are as follows. 

Two hundred twenty-five X-ray sources with a limiting luminosity of ~ 10 29 erg s _1 
are detected (§2J|, 202 of which are unambiguously identified with 2 MASS or UKIRT NIR 
sources ( §3.ip . These along with 4 other sources (identified as T-Tauri stars by their variable 
X-ray emission) are classified as stellar sources. The 19 remaining sources are likely AGN 



On the NIR color-magnitude diagram most of the detected X-ray stars occupy the locus 
of low-mass PMS stars. On the NIR color-color diagram 45 i^ s -band excess and 120 non-K s - 
band excess cluster members occupy the region with 5 mag< Ay <40 mag, but 10 stars are 
detected with low Ay and are identified as foreground field stars ( §3.2|) . The inferred X-ray 
properties are found to be consistent with IR properties of the stars through agreement with 
known L^-mass and Nu-Ay relations (55]). 

In the assumption of a universal XLF the XLF gives a lower limit distance to W40 
of 600 pc ( §6.ip . The XLF analysis shows that the total cluster population of PMS stars 
is >600 PMS stars assuming a distance of 600 pc, and a i^ s -band excess disk fraction of 
50% suggests a cluster age <1 Myr ( §6.2p . The IMF analysis is consistent with the XLF 
analysis on distance, total population, and disk fraction (§[7]). The X-ray observation is not 
sensitive to a protostellar population that might exist in W40 (jBontemps et al. 2010 ); only 
two protostellar X-ray candidates are identified ( §3.ip . 

The young age of the W40 cluster, inferred from the high disk fraction and the existence 
of gas in the cluster, and the ongoing star formation in the dust lane ( §9. 5[) indicate that the 
cluster should not be dynamically relaxed yet ( §9.ip . However, the cluster does show some 
signs of relaxation, including a nearly spherical morphology with no elongation ( §9.2.3p . a 
surface density profile that may be described by an isothermal sphere ( §9.2. ip . and mass 
segregation ( §9.3p . The cluster shape and profile may indicate the beginning of relaxation, 
but may also be geometric projection effects. Mass segregation d own to 1.5 M M in W40 (£19. 3p 



is unusual, but it is predicted in simulation studies of YSCs by iMoeckel fc Batd ( 120101 ) 



The W40 star cluster may be composed of populations with different ages ( §9.4p . Many 
of the stars without i^ s -band excess may belong to an older population which has progressed 
further toward gravitational relaxation, while many of the stars with i^-band excess, which 
are almost all located on the eastern side of the cluster, may be younger and may have not 
had time to mix with the rest of the cluster. It is likely that most of the high mass stars, 6 
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out of 8 having K s -band excess, are members of the second population, suggesting delayed 
massive star formation f£16.2D. 



The dust lane, which passes through the core and is detected by X-ray, optical, and 
IR absorption, may be the site of ongoing star formation ( §9.5p . In contrast, no evidence is 
found for star formation in the molecular core detected in 12 CO and 13 CO. The dust lane and 
molecular core have similar masses and may both be part of a ring of interstellar material 
surrounding the star cluster. However, the dust lane is undetected in CO, and might be an 
example of the destruction of molecular gas by stars or molecular gas freezing out on to dust 
grains. 

Diffuse X-ray emission is detected in W40 ( §10.2j) . However, it may be explained as 
X-rays from the unresolved stellar population. A powerful flare is observed from a ~4 M 
star without i^-band excess ( §10. 3p . 

A minority of studies consider W40 to be part of the Aquila Rift at a distance of 260 pc, 
rather than the more commonly adopted distance of 600 pc ( §l.ip . The lower distance would 
modify results in the following ways: the shape of the W40 XLF would be inconsistent with 
the COU P XLF; for W40 stars without K s -excess the shape of the IMF would be inconsistent 



with the (IChabrierl 120031 ) or COUP IMFs; the lower limit on the total cluster population 
would be 200 stars; only 2 stars would have masses above 10 M ; and mass segregation 
would be seen down to 0.3 M . 
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A. Photometric Mass Uncertainty 



We make an attempt to evaluate systematic uncertainty on masses using a testbed sam- 
ple of stars from the Taurus molecular cloud with known masses. The combined systematic 
and statistical error on masses of W40 stars is reported in Table [2j We do not attempt 
to estimate systematic uncertainty on stars with masses >5 M Q , and uncertainty on stars 
with masses up to 5 M is assumed to be similar to that derived from the Taurus sample 
with masses between 0.5 and 1.5 M Q . Sources of error on stellar mass include statistical 
photometric error, uncertainty on distance, uncertainty on median cluster age, intrinsic age 
spread, reddening law and magnitude of extinction, model uncert ainty including unmodeled 



effects of rotation, magnetic fields, and early accretion histor y (IBaraffe fc Chabrierl 12010 



Chabrier et al.l l2007t iMohanty et al.l l2009t iJeffries et al.ll2008l ). method uncertainty, bina- 
rity, and variability. 



A variety of evolutionary models are investigated by iHillenbrand fc White! (12004 ) who 
compare masses determined by placement on the Hertzspr u ng- Ru ssell (H-R) diag r am to 



dynamical masses for 27 PMS stars. For the paraffe et al.l ( 



998) and ISiess et all (119971 ) 



models used in our synthetic isochrones IHillenbrand fc White! ( 12004! ) found H-R masses are 
underestimated by ~10%. 

Ten out of the 27 PMS stars (7 single stars and 3 binaries) are from the Taurus star 
form ing region with a mean distance of 140 pc and a median age of 2.5 Myr (e.g. iGiidel et al. 
20071 ). The differences between dynamical mass and JHK S mass derived for these stars using 
the synthetic isochrone from 2.5 Myr and a distance of 140 pc, is used here to evaluate the 
magnitude of error on JHK S masses for W40. For binary stars we compare the JHK S 
mass to the mass of the most massive component. Table [S] provides dynamical masses for 
these stars, 2MASS photometry, inferred JHK S masses, and the percent error on the JHK S 
masses. 

The inferred difference in masses is around 30%, adopted here as an average systematic 
error on mass for W40 stars. However, this is likely a lower limit on uncertainty. The 
uncertainties from intrinsic age spread, model uncertainty, method uncertainty, binarity, 
and variability are partially accounted for assuming similarity with the testbed of sample 
stars from Taurus. Distance, median age, and reddening law and extinction effects on mass 
uncertainty in W40 are not accounted for. Distance is fixed at 600 pc and 260 pc (masses 
using 260 pc are a factor of 3-4 lower then masses using 600 pc), age is fixed at 1 Myr (masses 
would be 60% smaller usi ng 0.5 Myr and a factor o f 2 larger using 3 Myr), and the reddening 
law is assumed to be the iRieke fc Lebofskyl (119851 ) law. 
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Table 1. Basic Source Properties 





Source 




Position 










Extraction 










Charactcrist 


ics 




No 


CXOW40 


a (J2000.0) 


.5 (J2000.0) 


Error 


e 


NCt 


crNCt 


BGt 


NC h 


PSF Frac. 


Signif. 


Anom 


Var. 


Eff. Exp. 


ME 






(dag) 


(deg) 


(arcscc) 


(arcmin) 


(counts) 


(counts) 


(counts) 


(counts) 










(ks) 


(keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


160 


183130.20-020718.0 


277.875854 


-2.121672 


0.1 


2.1 


17.8 


4.8 


0.2 


17.9 


0.91 


3.3 





2 


34.5 


4.3 


161 


183130.50-020515.8 


277.877075 


-2.087735 


0.1 


1.8 


29.7 


6.0 


0.3 


23.8 


0.90 


4.5 





4 


34.6 


3.2 


162 


183130.56-020530.6 


277.877350 


-2.091831 


0.1 


1.7 


10.8 


3.8 


0.2 


10.9 


0.90 


2.4 





2 


33.0 


3.2 


163 


183130.73-020452.2 


277.878052 


-2.081161 


0.2 


2.0 


4.11 


2.8 


0.1 


0.9 


0.70 


1.4 





2 


34.6 


1.8 


164 


183130.77-015927.8 


277.878235 


-1.991042 


0.4 


6.7 


36.5 


6.8 


2.5 


28.4 


0.89 


5.0 





4 


27.1 


2.8 


165 


183130.84-020540.6 


277.878510 


-2.094621 


0.2 


1.8 


5.8 


3.0 


0.2 


4.9 


0.90 


1.6 





3 


33.3 


2.8 


166 


183131.01-020440.5 


277.879211 


-2.077922 


0.1 


2.2 


30.7 


6.1 


0.3 


24.8 


0.89 


4.6 





2 


34.5 


2.8 
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183131.69-020611.6 


277.882019 


-2.103231 


0.2 


2.0 


8.8 


3.5 


0.2 


4.8 


0.91 


2.1 
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34.4 


2.2 : 


168 


183131.89-020011.6 


277.882874 


-2.001029 


0.6 


6.2 


10.7 


4.1 


2.3 


10.5 


0.89 


2.3 





2 


31.0 


3.7 


169 


183131.89-020611.6 


277.882874 


-2.100822 


0.2 


2.0 


5.8 


3.0 


0.2 


4.8 


0.91 


1.6 





2 


34.4 


3.7 



CO 



Note. — Column 1: X-ray catalog sequence number, sorted by R.A. Column 2: IAU designation. Columns 3-4: Right ascension and declination for epoch (J2000.0). Column 5: Estimated 
standard deviation of the random component of the position error, \j a< ^. + a y ■ The single-axis position errors, a x and a y , arc estimated from the single-axis standard deviations of the PSF inside 
the extraction region and the number of counts extracted. Column 6: Off-axis angle. Columns 7-8: Net counts extracted in the total energy band (0.5—8 keV); average of the upper and lower la 
errors on net counts. Column 9: Background counts expected in the source extraction region (total band). Column 10: Net counts extracted in the hard energy band (2—8 keV). Column 11: Fraction 
of the PSF (at 1.497 keV) enclosed within the extraction region. A reduced PSF fraction (significantly below 90%) may indicate that the source is in a crowded region. Column 12: Photometric 
significance computed as net counts divided by the upper error on net counts. Column 13: Source anomalies: "1" — fractional time that source was on a detector (FRACEXPO from mkarf) is <0.9. 
Column 14: Variability characterization based on K-S statistic (total band): "2" — no evidence for variability (0.05 < Pks)'i "3" — possibly variable (0.005 < PkS *^ 0.05); "4" — definitely 
variable (-PxS < 0.005). Value "1" is reported for sources with fewer than four counts or for sources in chip gaps or on field edges. Column 15: Effective exposure time: approximate time the 
source would have to be observed at the aimpoint of the ACIS-I detector in Cycle 8 to obtain the reported number of net counts. Column 16: Background-corrected median photon energy (total 
band). Tablc^is published in its entirety in an electronic form. A portion is shown here for guidance regarding its form and content. 



Table 2. IR and Radio Counterparts to Chandra Sources 



No 


CXOW40 


2MASS 


Offset 


J 


H 


K s 


PhCCFlg 


MemFlg 


Av 


Mass 


UKIRT .Fig 


VLAJD 






( 


arcscc) 


(mag) 


(mag) 


(mag) 






(mag) 


M Q 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


160 


183130.20-020718.0 


18313020-0207180 


().() 


16.61± 0.15 


13.59± 0.04 


11.66± 0.03 


BAAccO 


W40-Ks 


20.09 


1.44±0.53 


1 


20 


161 


183130.50-020515.8 


18313049-0205158 


().() 


15.70± 0.05 


13.05± 0.04 


11.62± 0.03 


AAA000 


W40-Ks 


17.25 


1.48±0.47 


1 




162 


183130.56-020530.6 


18313056-0205306 


0.1 


15.57± 0.07 


13.09± 0.02 


11.82± 0.02 


AAAOOs 


W40-noKs 


15.68 


1.40±0.47 


1 




163 


183130.73-020452.2 


18313069-0204522 


0.5 


16.37± 0.10 


13.86± 0.04 


12.53± 0.03 


BAA000 


W40-noKs 


14.09 


0.58±0.19 


1 
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183130.77-015927.8 


18313077-0159275 


0.2 


16.72± 0.15 


14.34± 0.04 


13.24± 0.04 


CAA000 


W40-noKs 


13.38 


0.36±0.13 


L> 
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183130.84-020540.6 


18313084-0205409 


0.3 


16.04± 0.09 


13.51± 0.04 


12.30± 0.03 


AAA000 


W40-noKs 


14.72 


0.82±0.29 


1 
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183131.01-020440.5 


18313101-0204405 


0.1 


14.55± 0.03 


11.57± 0.02 


9.80± 0.02 


AAA000 


W40-Ks 


26.31 


12.27±0.35 


1 




167 


183131.69-020611.6 


18313168-0206114 


0.2 


17.03± • ■ ■ 


14.20± 0.04 


12.52± 0.03 


UAA000 


W40 


16.76 




1 




168 


183131.89-020011.6 














EG 











169 


183131.89-020611.6 


18313187-0206022 


0.7 


18.37± • ■ ■ 


14.78± 0.05 


13.44± 0.04 


UAA000 


W40 


25.37 




1 




Note. 


— Columns 1-2: X 


-ray source number and IAU designation. Column 


3: 2MASS source. Column 4: 


Chandra-2MASS positional 


offset. 


Columns 5-7: 


2MASS JHK S 


magnitudes. 


Column 


8: 2MASS photometry quality and confusion 


flag. Column 9: Membership flag: W40-noKS - W40 stars without K s 


-excess; W40-KS - W40 stars with if s -excess; FRG 


- foreground 


stars; EG - cxtragalactic source. Columns 10-11: Visual absorpt 


ion and stellar 


mass derived fr 


om the color-ma 


gnitude J vs 


J — H diagram assuming distance of 600 pc and age of 1 Myr. 


Column 


12: UKIRT flag: - 


no source; 1 - single source present; 1 


I - outside Field of View; 3 - double source present. Column 


13: VLA source (Rodr 


gucz. Rodney, 


&t Reipurth in preparation). 


Table [3] 


is published in its entirety in an electronic for 


m. A port 


on is shown here for guidance 


regarding its form and content. 
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Table 3. Additional Non X-ray Cluster Members 
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K s 




PhCCFlg 


A v 




Mass 






(mag) 


(mag) 


(maf 






(mag) 




M 




(1) 


(2) 


(3) 


(4) 




(5) 


(6) 




(7) 




18304958-0211176 


16.07±0.11 


14.62±0.06 


13.67±0.05 


BAA000 


6.43 


0. 


12 


±0 


.08 


18310021-0203383 


16.50±0.15 


14.16±0.04 


12.53±0.03 


CAA000 


12.67 


0. 


38 


±0 


04 


18311221-0206246 


16.46±0.13 


14.69±0.06 


13.54±0.05 


CAA000 


8.94 


0. 


15 


±0 


.05 


18311237-0212503 


16.48±0.14 


15.27±0.11 


14.34±0.09 


CBAOcc 


4.62 


0. 


07 


±0 


14 


18312171-0206416 


16.53±0.18 


14.53±0.10 


12.82±0.08 


CAA000 


10.30 


0. 


19 


±0 


.05 


18312271-0205484 


15.85±0.09 


14.42±0.05 


13.49±0.05 


AAAccc 


6.11 


0. 


13 


±0 


.06 


18312276-0204269 


16.01±0.09 


14.14±0.02 


12.97±0.02 


AAAcOO 


9.28 


0. 


23 


±0 


04 


18312330-0205139 


16.32±0.11 


14.04±0.07 


12.45±0.03 


BAAccc 


12.06 


0. 


38 


±0 


.39 


18312391-0201407 


15.89±0.08 


13.49±0.03 


12.00±0.02 


AAA000 


13.23 


0. 


68 


±0 


04 


18312396-0207040 


16.84±0.18 


14.93±0.07 


13.66±0.08 


CAAccc 


10.16 


0. 


15 


±0 


04 


18312403-0203282 


16.07± 0.08 


14.27± 0.04 


13.05± 


0.04 


AAA000 


8.79 


0. 


19 


±0 


.34 


18312554-0200332 


17.01± 0.18 


15.43± 0.10 


14.32± 


0.08 


CAA000 


7.69 


0. 


09 


±0 


.05 


18312611-0209246 


15.63± 0.05 


13.46± 0.03 


12.16± 


0.03 


AAA000 


11.19 


0. 


54 


±0 


22 


18312701-0212024 


16.51± 0.13 


14.06± 0.04 


12.53± 


0.03 


BAA000 


13.49 


0. 


46 


= 


.08 


18312701-0204530 


16.45± 0.12 


14.74± 0.07 


13.46± 


0.05 


BAAccc 


8.38 


0. 


14 


±0 


.04 


18312793-0204290 


16.14± 0.10 


14.42± 0.05 


13.36± 


0.04 


AAAccc 


8.38 


0. 


16 


±0 


.03 


18312846-0204187 


16.55± 0.13 


14.90± 0.06 


13.67± 


0.06 


BAA000 


7.95 


0. 


12 


= 


.05 


18312862-0209275 


15.68± 0.07 


14.18± 0.03 


13.22± 


0.04 


AAA000 


6.57 


0. 


16 


±0 


.03 


18312894-0209565 


14.80± 0.03 


12.80± 0.02 


11.61± 


0.02 


AAA000 


9.93 


0. 


74 


±0 


.06 


18312902-0206314 


16.94± 0.20 


14.98± 0.06 


13.70± 


0.07 


CAA000 


10.52 


0. 


15 


±0 


Hi 


18312919-0204551 


15.31± 0.04 


12.97± 0.02 


11.56± 


0.04 


AAAccc 


13.42 


1. 


05 


= 


.03 


18312938-0206490 


15.36± 0.05 


13.17± 0.04 


11.86± 


0.03 


AAA000 


11.37 


0. 


68 


±0 


.17 


18312972-0206022 


16.00± 0.09 


13.38± 0.04 


11.66± 


0.03 


AAA000 


15.62 


1. 


03 


= 


.05 


18313065-0208398 


16.45± 0.12 


14.92± 0.06 


13.76± 


0.04 


BAAccc 


7.00 


0. 


11 


±0 


.06 


18313331-0213012 


16.18± 0.10 


14.34± 0.04 


13.02± 


0.03 


AAA000 


9.11 


0. 


19 


±0 


.02 


18313457-0158555 


15.90± 0.08 


13.91± 0.04 


12.72± 


0.03 


AAA000 


10.00 


0. 


32 


±0 


.04 


18313566-0212534 


16.52± 0.11 


15.11± 0.08 


14.09± 


0.07 


BAAccc 


6.29 


0. 


09 


= 


.02 


18313595-0213305 


16.11± 0.09 


13.97± 0.04 


12.32± 


0.02 


AAA000 


10.90 


0. 


34 


±0 


11 


18313646-0209212 


14.37± 0.04 


12.30± 0.04 


10.97± 


0.02 


AAA000 


12.09 


1. 


40 


±0 


.03 


18313789-0209080 


16.75± 0.13 


15.39± 0.10 


14.35± 


0.12 


CABccc 


5.83 


0. 


08 


±0 


.03 


18313791-0207453 


16.34± 0.10 


14.44± 0.04 


13.23± 


0.03 


BAA000 


9.75 


0. 


19 


±0 


.23 


18313815-0214236 


16.80± 0.14 


14.77± 0.06 


13.48± 


0.05 


CAA000 


10.91 


0. 


17 


±0 


.06 


18313832-0211449 


16.05± 0.07 


14.70± 0.04 


13.82± 


0.04 


AAA000 


5.62 


0. 


11 


±0 


.03 


18313878-0211155 


16.45± 0.10 


14.85± 0.06 


13.76± 


0.06 


BAAOcO 


7.63 


0. 


12 


±0 


.03 


18313979-0203558 


16.99± 0.17 


15.29± 0.10 


14.10± 


0.07 


CAA000 


8.53 


0. 


10 


±0 


43 


18314106-0211475 


16.73± 0.15 


15.65± 0.12 


14.72± 


0.12 


CBBOOO 


3.58 


0. 


06 


±0 


14 


18314362-0209545 


16.54± 0.10 


15.38± 0.08 


14.28± 


0.07 


BAAOOO 


4.15 


0. 


07 


±0 


.23 


18315335-0205248 


17.02± 0.19 


15.43± 0.08 


14.34± 


0.08 


CAAOcc 


7.77 


0. 


09 


±0 


10 


OS 3a 


8.91± 0.03 


8.02± 0.05 


7.37± 


0.02 


AAAOdd 


8.84 


15.00 


±0.72 



CO 



Note. — Column 1: 2MASS source. Columns 2-4: 2MASS JHK S magnitudes. Column 5: 2MASS 
photometry quality and confusion flag. Columns 6-7: Visual absorption and stellar mass derived from the 
color-magnitude J vs. J — H diagram assuming distance of 600 pc and age of 1 Myr. 



Table 4. X-ray Spectroscopy 
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XPHOT 






No 


CXOW40 


JVC 


ME 


MemFlg 


l°g( L hl) 


log(itl) 




rog(JVffl) 




log(itcl) 


log(N H2 ) 


!og( L fic2) 


log(Ltc2) 






(cnts) 


(keV) 




(ergs s -1 ) 


(ergs s -1 ) 


(keV) 


(cm- 2 ) 


(ergs s -1 ) 


(ergs s -1 ) 




(cm- 2 ) 




(ergs s- 1 ) 


(ergs s 


l ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 




(12) 




(13) 


(14) 




160 


183130.20-020718.0 


17.85 


4.30 


W40-Ks 


29.79 


29.79 


1.73 


23.26 


30.52 


31.02 


23 


.15 ± 0. 


19 


30.31 ± 0.17 


30.73 ± 


24 


161 


183130. 50-020515. S 


29.74 


3.22 


W40-Ks 


29.86 


29.89 


5.03 


22.46 


29.99 


30.22 


22 


.64 ± 0. 


19 


30.07 ± 0.14 


30.52 ± 


.27 


162 


183130.56-020530.6 


10.77 


3.23 


W40-noKs 


29.47 


29.47 


1.53 


22.96 


29.99 


30.57 


22 


.71 ± 0. 


26 


29.74 ± 0.20 


30.22 ± 


.36 


163 


183130.73-020452.2 


4.88 


1.76 


W40-noKs 


28.62 


28.86 


0.55 


22.54 


28.99 


30.65 














164 


183130. 77-015927.S 


36.54 


2.84 


W40-noKs 


29.88 


29.92 


1.23 


22.71 


30.27 


31.00 


22 


.50 ± 0. 


15 


30.09 ± 0.12 


30.56 ± 


.25 


168 


183130.84-020540.6 


5.82 


2.76 


W40-noKs 


29.09 


29.12 


1.72 


22.66 


29.38 


29.88 


22 


.52 ± 0. 


28 


29.25 ± 0.23 


29.76 ± 


.34 


166 


183131.01-020440.5 


30.73 


2.83 


W40-Ks 


29.75 


29.79 


1.76 


22.56 


30.00 


30.50 


22 


.50 ± 0. 


14 


29.99 ± 0.12 


30.41 ± 


24 


167 


183131.69-020611.6 


8.78 


2.17 


W40 


29.11 


29.18 


1.83 


22.37 


29.27 


29.74 


22 


.30 ± 0. 


26 


29.23 ± 0.21 


29.72 ± 


211 


168 


183131.89-020011.6 


10.67 


3.72 


EG 


29.48 


29.48 


1.39 


22.97 


30.04 


30.67 


22 


.88 ± 0. 


28 


29.88 ± 0.21 


30.33 ± 


.32 


169 


183131.89-020611.6 


5.78 


3.66 


W40 


29.13 


29.15 


1.83 


22.75 


29.46 


29.93 


22 


.89 ± 0. 


34 


29.57 ± 0.28 


30.10 ± 


40 



o 



Note. — Columns 1-5: For convenience, these columns reproduce the source identification, net counts, median energy, and membership flag from Tables |T| and [2l Columns 6-11: Apparent 
X-ray luminosities, plasma temperature, column density, and intrinsic X-ray luminosities inferred from XSPEC fitting. Columns 12-14: X-ray column density, intrinsic luminosities, and their 
errors (summed in quadrature statistical and systematic errors) obtained with XPHOT. pro. X-ray luminosities are derived assuming a distance of 600 pc. Tablc[^is published in its entirety 
in an electronic form. A portion is shown here for guidance regarding its form and content. 
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Table 5. Completeness Limits and Scale Factors 





Hard-Band XLF 




Total-Band XLF 




IMF 




Sub population 


log Completeness Limit 


Scaling 


log Completeness Limit 


Scaling 


Completeness Limit 


Scaling 




(erg s _1 ) 




(erg s _1 ) 




(M ) 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


all w40 


30.2 


0.70 


30.6 


0.70 


1.5 


0.72 


no /f s -band excess 




0.30 


30.0-30.4 


0.35 


0.5 


0.32 


i-C s -band excess 




0.30 


30.8 


0.35 


1.5 


0.40 



Note. — Column 1: Population examined. We examine the entire sample of W40 cluster members, cluster members with Jf s -band 
exc ess, and c luster members without it s -band excess. Columns 2,4,6: The completeness limits at which the scaled COUP XLF model 
and [Ch abricr I 2003) model depart from the W40 XLF and IMF, respectively. Columns 3,5,7: The scale between the model needed to 
fit the COUP sample of > 800 stars and the model needed to fit the W40 subpopulation. 
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Table 6. NIR and X-ray Properties of Early- Type Star Candidates 



No 


Name 


Arc 


Mass 


log(iBol) 


if s -band excess 


N H 


kT 


log(itc) 






(cnts) 


(M ) 


(L Q ) 




(10 22 cm" 2 ) 


(keV) 


(ergs s -1 ) 


(1) 


(2) 


(3) 


(4) 


5 


(6) 


(7) 


(8) 


0) 


46 


183117.12-020911.5 


9 


12 


4.11 


NoKs 


2.32 


1.64 


30.06 


99 


183124.00-020529.5, OS 2a 


58 


30 


5.25 


Ks 


3.99 


1.24 


30.95 


122 


183126.02-020517.0 


1177 


16 


4.46 


Ks 


1.69 


3.55 


31.63 


141 


183127.84-020523.5, OS la 


145 


<32 


<5.36 


Ks 


1.78 


1.21 


30.94 


144 


183128.01-020517.1 


11 


14 


4.24 


Kb 


1.19 


0.68 


29.84 


145 


183128.67-020529.8 


417 


21 


4.85 


Kb 


3.41 


3.70 


31.42 


166 


183131.01-020440.5 


41 


12 


4.09 


Ks 


9.38 


0.90 


31.47 




OS 3a 




15 


4.36 


NoKs 









Note. — Columns 1-2: X-ray source number and IAU designation. Column 3: X-ray source net counts. Columns 4-5: NIR 
photometric stellar mass and bolometric luminosity. Column 6: K s -band excess indicator. Columnns 7-9: Inferred X-ray spectral 
properties. 
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Table 7. Cluster Structure 



Subpopulation 

(i) 


Center Distance 
(pc) 
(2) 


Center Angle 
(degrees) 
(3) 


Median Radius 
(PC) 
(4) 


ALL 


0.00 


NA 


0.50 


LM 


0.02 


319 


0.60 


IM 


0.01 


270 


0.36 


HM 


0.12 


107 


0.08 


NKE 


0.11 


273 


0.59 


KE 


0.12 


80 


0.33 


LM; NKE 


0.13 


276 


0.66 


IM; NKE 


0.08 


232 


0.41 


LM; KE 


0.27 


62 


0.54 


IM; KE 


0.04 


348 


0.35 



Note. — Column 1: The subpopulation of stars used for structure anal- 
ysis; LM - low mass stars (_M < 1.5 Mq); IM - intermediate mass stars 
(1.5 < M < 10 M Q ): HM - high mass stars (M > 10 Mq); NKE - 
stars with no _?C s -band excess; KE - stars with i^-band excess. Column 
2: Distance from the cluster center to the subpopulation center (median 
position). Column 3: Position angle of the subpopulation center with re- 
spect to the cluster center. Column 4: The radius of the circle around the 
subpopulation center that contains half the stars in the subpopulation. 
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Table 8. Dynamical and JHK S Masses of Taurus Stars 



Name 


J 


H 


K B 




MjHK B 


% Error 




(mag) 


(mag) 


(mag) 


(M ) 


(M ) 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


MWC 480 


5.41 


5.34 


5.34 


1.65 a 


2.68 


62 


BP Tau 


9.10 


8.22 


7.74 


1.32° 


1.22 


-7 


LkCal5 


9.42 


8.60 


8.16 


0.97" 


0.73 


-24 


GM Aur 


9.34 


8.60 


8.28 


0.84" 


0.66 


-21 


DL Tau 


9.63 


8.68 


7.96 


0.72 a 


0.78 


8 


DM Tau 


10.44 


9.76 


9.52 


0.55" 


0.30 


-44 


CY Tau 


9.83 


8.97 


8.60 


0.55" 


0.58 


6 


045251+30168A 


8.96 


8.32 


8.13 


1.45 c 


0.98 


-32 


045251+30168B 


8.96 


8.32 


8.13 


0.81 c 






UZ Tau Aa 


8.45 


7.60 


7.35 


1.02" 


2.10 


105 


UZ Tau Ab 


8.45 


7.60 


7.35 


0.294° 






GG Tau Aa 


8.67 


7.82 


7.36 


0.67 ad 


1.87 


180 


GG Tau Ab 


8.67 


7.82 


7.36 


0.61 ad 






Note. — Column 


1: Star 


name. Col 


umns 2-4: 


2MASS 


J H K s mag 


nitudes. Col- 



umn 5: Dynamical masses. Column 6: 2M ASS photometric mass es. Column 7 : 
P ercent error in phot ome tric masses. a ISimon et al . 12000). '' D ut re v et al ] i 20031 ). 
c ISteffen et al.l l200l]) . d IWhite et al.l il999l) . 
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I (deg) 

Fig. 1. — (a) 5°x5° DSS Red image of the neighborhood of W40, including Serpens South 
and Serpens Main. Nearby dark clouds and star forming regions are indicated by circles. 
The Chandra ACIS-I field of view is indicated by a white square, (b) Enlarged image of the 
ACIS-I field of view. Nebulous optical emission and a few stars projected on W40 may be 
seen. 
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Fig. 2.— (a) 24' x 24' 2MASS i^-band and (b) 11' x 11' UKIRT K-band images of the 
W40 cluster. The solid line indicates the Chandra ACIS-I field of view and the dashed line 
indicates the UKIRT field of view, (c) Adaptively smoothed 17' x 17' ACIS-I image of the 
W40 cluster. Smoothing has been performed on the total 0.5-8.0 keV band at the 2.5 a level, 
and the grey scale is logarithmic. Most of the 225 X-ray sources are visible, (d) The 6' x 5.6' 
center of the previous image is shown in greater detail. 
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Fig. 3. — NIR color-color (a) and color- magnitude (b) diagrams for 164 X-ray sources clas- 
sified as cluster members with K s (triangle) and without -?T s -band excess (circle), and fore- 
ground stars (square). In panel (a): the solid curved line shows sites of intrinsic JHK S colors 
of 1 M yr (see ^3.21) stars ; the s traight solid line is a reddening vector (using the reddening 
law of iRieke fc Lebofskyl (119851 )) originating at M — 0.2 M ; the dashed line is the CTT 



locus ( Meyer et al. 



19971 ). In panel (b) the solid line is 1 Myr isochrone, and the dashed lines 



are reddening vectors of Ay = 10 mag shown for various star masses. 
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Fig. 4. — (a) X-ray luminosity as a function of mass for lightly obscured PMS population 
in COUP (grey x), and W40 stars with (triangle) and without (circle) /A^-band excess, 
(b) X-ray absorbing column density from the X-ray spectral analysis with XPHOT and 
visual absorption fr om the IR p hotometry. The curve gives the gas-to-dust relationship 
N H = 2.2 x 10 21 A V fevterlll99( 
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Fig. 5. — Comparison between XLFs of the ONC and W40 populat ions. The upper blac k 
histograms shows the COUP unobscured cool sample of > 800 stars (jFeigelson et al.ll2005l ): 
the red and blue histograms show samples of W40 stars with and without ^-excesses, 
respectively. Error bars indicate 68% confidence intervals (1 a) from Monte Carlo simulated 
distributions when X-ray luminosities are randomly drawn from Gaussian distributions with 
the source's measured logLx values and estimated errors AlogLx- The green lines are 
added to aid the eye; they are based on the shape of the ONC XLF and scaled downward to 
match the W40 population without _fT s -excess. Hard-band and total-band XLFs are shown 
on the left and right, respectively. Panels (a)-(b), (c)-(d), and (e)-(f) are for the assumed 
distances to W40 of 260, 600, and 900 pc, respectively. 
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Fig. 6. — Comparison between IMFs of the COUP (dashed lines) and the W40 populations 
(solid lines). The top two panels show IMFs for 600 pc for stars without K s -band excess 
(a) and with i^ s -band excess (b). The bottom two panels show IMFs for 260 pc for stars 
without K s -band excess (c) and with i^-band excess (d). Bin sizes are adjusted to cover 
the entire region of JHK S photo metric mass deg eneracy with a single bin ( §3.21) . and 1 a 
simulated errors are shown. The IChabrierl (120031 ) IMF (dotted lines), integrated to match 
the binning, is scaled to the COUP and W40 IMFs. 
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Fig. 7. — Total-band (0.5-8.0 keV) X-ray luminosity as a function of bolometric luminosity. 
Massive star candidates are circles and lower mass stars are triangles. Dotted lines indicate 
L x /L Bol = 10-\ 10- 2 , 10" 3 , 10~ 4 , 10" 5 , 10" 6 , 10" 7 , 10" 8 , and 10~ 9 . 
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Fig. 8. — The surface density of W40 cluster members detected in the X-ray, smoothed with 
a 1.5 arcmin Gaussian kernel. Contours are drawn starting at 100 stars degree -2 and at 
every subsequent increase by a factor of 1.5 in surface density. 
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Fig. 9. — (a) The radial profile of all W40 cluster members (black squares), stars with K s - 
band excess (blue squares), and stars with no i^ s -band excess (red squares). y/N error bars 
are used. On top, the radial profile of stars in COUP is plotted for comparison (black solid 
line), and it is scaled to match the various W40 radial profiles, (b) The radial profile of all 
W40 cluster members is shown in black, and the best King profile ( King Ill962l ) found from 
fitting the cumulative distribution is shown in blue. 
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Fig. 10. — The minimal spanning tree of W40. Chandra-detected cluster members arc 
vertices. 
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Fig. 11. — (a) The cumulative distribution function of azimuthal angles, (b) The solid line 
marks the median distance of stars from the cluster center in each 60 degree segment. Error 
bars are plotted, and the best fitting constant radius is shown (dashed line). The boundary 
of the center of the cluster from which stars are excluded is shown with a dotted line. Stars 
are indicated by diamonds. 
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Fig. 12. — Positions of W40 stars with (black triangles) and without (grey circles) K s - 
excesses. The 17' x 17' ACIS-I Fiel d of View (grey box) is shown. The contour marks the 
13 CO emission from lZhu et al.l (120061 ). Panels show different stellar mass strata: M < 1.5 M 
(a), 1.5 < M < 10 M Q (b), and M > 10 M (c). 
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Fig. 13. — (a) Positions of Chandra-detected K s -excess stars (triangles) and non-fT s -excess 
stars (circles) and addition al iC-excess star s not detected in by Chandra (squares). The 
13 CO contour is also shown ( jZhu et al.l 120061 ) . (b) The cumulative distributions of Chandra- 
detected stars with and without fC,-excess are shown. 
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Fig. 14. — (a) Positions of the W40 stars with J — H > 2.2 mag (blac k circles) and J — H < 
2.2 mag (grey circles). The contour marks the 13 CO emission from IZhu et al.l ( 120061 ). (b) 
Positions of the W40 stars with J — H > 2.2 mag over-plotted on the red Digitized Sky 
Survey image with inverted colors. The box outlines the Chandra ACIS-I field of view. 



Fig. 15. — The smoothed Ay image is shown with intensity ranging from log(Ay) =0.7 to 
1.3 mag. Star positions are indicated by black triangles for stars with i^-band excess and 
grey circles for stars wit hout ^-band e xcess. Contours show the molecular core in 13 CO 
J = 2 — 1 emission lines (jZhu et al.ll2006l ). A white polygon outlines the region with greatest 

Ay. 



60 




Fig. 16. — (a) Total-band (0.5 — 8 keV) Chandra ACIS-I image after removal of point 
sources. This image is adaptively smoothed with a Gaussian kernel and normalized by the 
smoothed exposure map. The central polygon marks the extraction region for the diffuse 
X-ray emission with the carved out central part due to the contribution of the PSF wings of 
the point sources in the core of the cluster. Three background regions for the diffuse emission 
are also indicated, (b) Comparison of the composite spectra of the W40 stars, AGN, and 
the diffuse emission. 



